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AN OCEAN WAVE MEASURING BUOY 


ABSTRACT 

We describe the design, construction, and performance of 
a surface wave-following buoy that can measure the directional 
spectrum of 0.06 to 0.5 Hz ocean waves with an accuracy of 10%, 
and an angular resolution of around 90°. The buoy is uisc 
shapeo, 1.5 m (5 feet) in diameter, weighs 150 kgm (300 lb), 
and is completely self-contained. Wave height is measured by 
an accelerometer mounted on the inner gimbal of a vertical gyro, 
wave slopes by the tilt of the buoy about the gyro's vertical 
axis, ana buoy heading by a gyro stabilized compass. Data from 
the transducers are digitized and recorded on computer compa¬ 
tible magnetic tape. Wave spectra are calculated from the data 
by computer. The accuracy of the measurements is verified by 
calibration in a wave tank, and by the internal consistency of 
the data. 

1. INTRODUCTION 

The wave measuring buoy described in this report was 
designed to provide directional spectra ">f ocean waves for 
comparison with radar scatr.er data. Initially we sougnt to 
purchase a buoy, but a quick survey of available instruments 
showed that none were capable of measuring a directional 
spectrum. Thus we were forced to design and construct one. 
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A number of techniques exist for determining the direc¬ 
tional distribution of a wave field. Essentially, the measure¬ 
ment requires a conerent sample of the wave field over many 
wavelengths. The directional resolution is proportional to tne 
number of wavelengths used. The measureme it of the ocean wave 
field is particularly difficult because it is difficult to 
establish a fixea reference point to //hich wave heights can be 
related. Wave data can be measured at several places using an 
array on a stable platform, or, alternately, the surface eleva¬ 
tion and two components of slope (tilt) can be measured at one 
point using an inertial reference. 

A buoy that measures wave height and tilt (commonly 
called a pitcn-and-roll buoy) is particularly simple, ard can 
be small and easily handled. We have chosen this technique. 

Several pitch-and-roxl buoys have been built in the past. 
The first were designed and built by the National Institute of 
Oceanography in England. One was contained in a 5 feet 6 inch 
cast aluminum, ellipsoidal hull (Longuett-Higgins, Cartwright, 
and Smith, 1963). Later versions used 54" torus shaped hulls. 
Still later, another version was built by Hudson Laboratories 
of Columbia University (Saenger 1969a, b; Goldberg anil Goldberg, 
19t>9; Jordan, 1969). All of these buoys required a ship to 
stand by for recording data and supplying power. Our design 
borrows heavily from these previous instruments and is an 
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extension of their concept: It is completely self contained 
and can operate unattended. 

2. PRINCIPLE OF OPERATION 

The hull of a pitch-and-roll wave measuring buoy is disc 
shaped, of shallow draft, and is radially -.ymmetric about an 
axis normal to the water surface. Such a hull will follow the 
surface of a wave provided the wavelength is sufficiently long; 
and its high frequency response will be independent of wave 
direction. The buoy contains an accelerometer mounted on the 
axis of a vertical gyro and a gyro stabilized compass. The 
vertical component of acceleration, when doubly integrated, 
gives the sea surface height. The buoy tilts, when referenced 
to the vertical axis of the gyro and to the compass heading, 
gives the wave slopes in a fixed (North certered) coordinate 
system. These variables can be related to the ocean-wave 
directional spectrum and its lower-order moments. 

Let 


U-L, c 2 , c 3 ) = U, H/3x, a;/3y) (2.1) 

be the wave height and slopes measured by the buoy as a func¬ 
tion of time in a coordinate system with x,y pointing North and 
East (Fig. lb). The spectrum F^(k, ,k ) of c^(x,y,t Q ) at some 
instant t Q is so defined that 
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Similarly, for time series ;^(x,y;t) at one point x o y D , 
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where u> is a radian frequency, k is the wave number. 

By the ergodic theorem the space and time averages can be 
equated: 


F , (ui, 8) dui = F. (k, B) k dk . (2.4) 

l i 

For deep-water waves, w 2 = gk, wnere g is the acceleration 
of gravity, and the required Jacobian is 


F i (k,s) 

F^ (u), 6) 


i <L_ = 

k dk 2w 3 


(2.5) 


It is convenient to express all measured spectra ir terms 
of F. (k ,K ), F. (ui, 0) , etc., e.g. the contributions (per unit 

■L A y X 

wavenumber space, per unit frequency-radian, etc.) to the mean- 
square surface elevation. We omit the subscript "i" and refer 
to F ( ) as simply the wave spectrum. 

The Cartesian moments are written 
00 
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where 
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For an elementary wave train, the wave height and slope are 
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The co-spoctra C i j and quadrature-spectra j of any pair of 
quantities ; i and can be expressed in terms of the moments 
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Furthermore, a trigonometric identity gives 
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Equations (10) and (11) serve 
buoy data. 


to estimate the accuracy of the 
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The five moments determined by the buoy can serve to 
evaluate the first five Fourier terms of the directional 
distribution of ocean-wave energy (at each frequency) 

F< 5 )(k,6) = 1/2 a Q + a^cosB + b^inB + 

a2COs2B + b2si. 26 


with 


(a Q ,a^,b^,a 2 /b 2 ) 


(N 00' N 10' ,) Ol' N 20- N 02' 2N il ) !2 ' 12) 


The terminated Fourier expansion can be expressed in the 
form (Longuett-Higgins, Cartvright/ and Smith/ 19 63) 

F (5) (k,B) = i_J 2 F(k,B)W(B'-B) dB' (2.13) 

2 it o 

where W = 1 + 2cos(6'-6) + 2cos2(6'-6) can do regarded as a 
weighting function associated with the buoy measurements. When 
B'-B * 0, +44°, W = 5, 5/2 respectively. Thus the angular 
resolution of the tilt buoy can be taken at 88°. This is net 
very good, but it is the penalty one must pay for having a 
simple instrument. 


3. DESIGN PARAMETERS 

The wave buoy was to be used primarily in a trade-wind sea. 
To meet our requirements it must measure the ocean wave spectrum 
in the frequency band of 0.06 to 0.5 Hz (2-16 sec period) while 
wind speeds range up to 15 m/sec. Furthermore, it must be com¬ 
pletely self-contained and operate *ror up to 24 hours unattended. 


Typically, we expect to place it in operation and than conduct 
other experiments nearby. After tnree or four hours we would 
return and recover the buoy. The measured spectra should have 
an error of no more than 10% in amplitude or 10° in direction. 
These criteria determine the buoy size and transducer accuracy. 

To estimate the accuracy required of the transducers using 
these criteria we assume an ocean wave spectrum proposed by 
Pierson and Moskowitz (1964) 


S (w) = (ag 2 (d -5 ) exp [-Bg 4 (vid) " 4 ] 


(3.1) 


where 


2 IT 

3 (u) = / F((d,B) d B 

0 

and v is the mean wind velocity. We use a = 8.1 x 10 and 
B = 0.74. The acceleration spectrum derived from (3.1) is: 


S ((d) = ag 2 to -1 exp [-Bg 4 (V(d) “ 4 ] (3.2) 

A 

The root-mean-squc-'e acceleration is: 

o° 

<£ o 1 / 2 = (/ S,(<d) da,) 1 / 2 (3.3) 

A o A 

This integral diverges because of the contributions at large 
w. However, a buoy will respond only to waves whose frequency 
is less than some frequency, 0; its acceleration is found 
by integrating (2.3) to this upper limit. This gives: 




1/2 


(-cg 2 /4 E i (-Bg 4 (vO) -4 ] ) l / 2 


(3.4) 
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Here E ± is the logarithmic integral defined by Jahnke and Emae 
(1945, p. 1). We expect the buoy to respond to waves whose 
wavelength is roughly twice the buoy diameter or greater (see 
the end of this section). Waves at the high frequency cut-off 
typically have a frequency of one Hertz. Letting 1 = 2-/sec, 
>1/2 = 13% q when v = 15 m/sec and 5% g when v = 2 m/sec. 

To measure these accelerations with an accuracy of 10% requires 
an accelerometer with a total error band of less than 0.5% g 
over a range of 0-2 g. 

The directional spectrum (2.12) is calculated from the 
slope spectra and co-spectra. The accuracy is limited mainly 
by the coefficients a 2 and b 2 which depend on the two slope 
measurements. Each must be measured to an accuracy of 10% + /2 
if the two components, taken together, are to have a total 

error of less than 10% on average. 

The root-mean-square slope is directly related to the 
acceleration spectrum through the dispersion relation: 

<c 2 + ^ 2 > = k 2 <; 2 > = ^ 3 ’ 5) 

The acceleration measured in units of g is the slope in 
radians. A 2 m/sec wind gives an RMS surface slope of 3°, a 
15 m/sec wind gives 7°. A vertical gyro capable of measuring 
these angles to an accuracy of 7% must have an error of less 
than 0.2° over the band ±30°, the Stokes limit for a progres¬ 


sive wave of maximum amplitude. 
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The operation of the wave buoy requires a hull which 
accurately follows the water surface. One that is disc shaped 
and of shallow draft will perform well until the wavelength 
becomes less than about twice its own diameter. Since a 2 sec 
wave has a wavelength of 6.2 meters, the buoy should be some¬ 
what smaller than three meters. On the other hand, it should 
be easily handled. As a compromise, we arbitrarily chose a 
hull 1.5 m (5 ft) m diameter, with a draft of about 8 cm. 

Such a buoy displaces 150 kgm. Hudson Laborator.es' experience 
with smaller buoys indicated that a hull of this size will sur¬ 
vive 15 m/sec winds with only a small chance of being upset (of 
course it would not be sunk). 

We estimate the response of this size buoy hull from Kim 
(1966). From his figure 11, the half-power point in heave 
response occurs when a = Ru 2 g _1 = 1.5, R is the radius of the 
buoy (0.75 m). This occurs at a frequency of 0.7 Hz. In a 
similar manner, the half-power point in the pitch response 
occurs when a = 3.8 (Kim's figure 12). This corresponds to a 
frequency of 1.0 Hz. The hull responds to wave slope slightly 
better than to wave height near f a cut-off frequency. In 
either case, the response is within 10% of unity at frequencies 
of 0.5 Hz and lower, so the hull should meet our requirements. 

4. BUOY CONSTRUCTION 

The buoy consists of four main systems: 1) the trans¬ 
ducers, 2) the signal conditioning and recording system, 3) 
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safety and recovery system, and 4) the buoy hull that houses 
and supports this equipment. We will discuss each in order. 

The sources of supply and cost of equipment discussed below, 
together with detailed electrical schematics, are included in 
the appendix. 

The transducers consist of a vertical gyro, an acceler¬ 
ometer mounted on the inner gimbal of the gyro, and a gyro 
stabilized magnetic compass. 

The gyro-accelerometer assembly was made by Honeywell, 
and consists of a vertical gyro (part number JG7044A45) and a 
quartz-fiber accelerometer (GG326C1). The accelerometer has a 
range of 0 to 2 g, a total error band of 0.01% g, a drift in 
sensitivity of 0.02% g/°C, and n zero point of 4 x 10“ 6 g/°C. 
The gyro gimbal remains vertical to within +0.25°, and has a 
full scale range of ±30°. The tilt of the buoy is measured 
with a resistance potentiometer having a 1% linearity. The 
entire assembly is 26 x 15 x 14 cm in size, weighs 4 kgm, and 
uses about 50 watts of power. 

Buoy heading is measured by a Humphrey North Seeking Gyro 
(DG04-0122-1). This consists of a gyro stabilized and gimbaled 
magnet whose position is measured by a resistance potentiometer. 
The unit has an accuracy of ±1°, and is linear within ±1%. 
Physically, it is 21 x 8 x 8 cm, weighs 1 kgm, and uses 10 
watts of power. 


The recording system conditions thp signals from the 
transducers, converts them to digital numbers, and records 
them on computer compatible magnetic tape. In addition, it 
controls the operation of the buoy, times its operation, and 
writes the time on the tape. The system was designed and built 
by Monitor Laboratories. The tape recorder is a standard unit 
manufactured by Precision Instruments (PI1387) . The assembly 
is remarkably small (see figure 2), and its standby power 
requirement was less than two watts. 

The signals from the transducers vary between +5 volts, 
and contain some noise at higher frequencies, particularly at 
400 Hz (the frequency of the gyro supplies). To reduce alias¬ 
ing errors, the signals are sent through low-pass filters. The 
voltages from the potentiometers (tilts and heading) go through 
a simple RC filter with a cut-off frequency of 20 Hz. The 
accelerometer signal is more severly attenuated since it must 
be integrated twice to obtain sea surface heights. Even a 
small amount of aliased power could cause severe error:; in the 
measurement of low-frequency wave heights. Ihis signal goes 
through a two-pole filter with a cut-off frequency of 1.0 Hz. 
The complex filter response functions are: 

L = (w 2 - 1) (w 4 + l)' 1 + j/2 uU 4 + l)' 1 
= (w 2 + l) -1 + jw(w 2 + l) -1 


L 


(4.1) 












wnere uj is frequency in Hz. Measured responses are very close 
to values calculated from these equations: within 4% in ampli¬ 
tude and 1/2° in phase at the half-power frequency. 

The filtered signals are multiplexed, converted to 12 bit 
binary numbers (an accuracy of approximately ±0.01%), and 
recorded on 1/2" magnetic tape. The recording format i, vari¬ 
able and controlled by switches. The correct time (year, day, 
’-'nr, minute) is recorded at the beginning of each record. The 
time base was a quartz crystal having an accuracy of +4 seconds/ 
day over a 50°C temperature range. 

A 7 track, incremental tape recorder records the data at 
a density of 200 characters per inch and at a rate of up to 200 
steps per second. It is rugged, can operate with ? g RMS accel¬ 
erations, uses one watt of power during standby and 40 watts 
while recording. The total amount of energy required to write 
600 feet of tape is fixed at approximately 7 amp-hours from a 
12. volt supply (72 joules). The rate of energy used is deter¬ 
mined by the rate data is recorded. Typically, recording 144 
bits/second uses 4 watts (144 bits/second = four 12-bit words, 
three times a second). 

The recording system controls the buoy operation. The 
correct time is entered in the clock by switches. Once started, 
the system turns on power to the equipment 10 minutes before 
the next hour, and begins recording on the hour. Data are 
recorded in the format selected. After the selected amount has 
been written, recording stops and power is turned off. One, 
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two, or four hours later the cycle repeats, thus allowing the 
buoy to operate unattended for a number of hours. 

Typically, each transducer is sampled 3.125 times/sec; 2048 
samples are written in each record; and 16 records comprise a 
file. The records contain 10 minutes of data, the file about 3 
hrs. Each file uses about 9 meters (27 feet) of tape. Spectra 
calculated from this data have a resolution of 1.5 x 10“ 3 Hz and 
a Nyquist frequency of 1.6 Hz. This format ensures that low 
frequency waves are adequately resolved, aliasing errors are 
small, and the time series are easily handled in the computer. 

Power for the operation of the gyros and recording system 
comes from sealed, gel-cell, lead-acid batteries. In normal 
operation the buoy uses 7 amperes at 12 volts. Twenty-four 
hours of continuous operation requires 168 amp-hours of energy, 
but intermittent operation reduces this considerably. The buoy 
has six 20 amp-hour batteries. They eperate in any position, 
and need not be protected from sea water. Each battery is iso¬ 
lated by diodes: failure of one does not affect the operation 
of the others. The batteries weigh 48 kgm and are a major part 
of the buoy weight, but are small in size (0.02 in 3 ). 

The 12 VDC power is converted into other voltages required 
by the systems (see figure 3) . The gyros operate from 115 and 
26 VAC, 400 Hz. The first voltage is produced by a switched 
inverter (Nova). Switching transients occur on the 12 VDC 
lines and are reduced by the input and RF filters. The 26 VAC 
is derived from a transformer. A small amount of +15 and +5 




































/DC is produced by twc Acopian power supplies. This serves the 
potentiometers, accelerometers, and data recorder. 

The buoy is designed to operate unattended, consequents 
it has systems that help in finding it once it is out of sight. 
It also recognizes and responds to certain dangerous conditions 
such as overturning in large seas. In this event, the gyros 
tumble and no data can be obtained. Other mishaps include loss 
of power (batteries discharged), broken magnetic tape, or end 
of tape. The occurrence of any of these events turns off tne 
power to the gyros, terminates the recording of data, and causes 
the buoy to call for help on its radio. 

To aid in finding the buoy, it is painted bright yellow, 
has a flashing zenon light, and a citizen band (27 MHz) radic 
transmitter. The buoy can be seen fox several hundred meters 
in 2 meter seas. The flashing light can be seen for one kilo¬ 
meter at night. The radio can be heard from 10 km away. The 
transmitter sends out 500 m watt pulses and has two antennas, 
one on top and one on the bottom of the hull. In the event of 
a mishap the pulse rate doubles. If the buoy overturns, the 
bottom antenna is activated. The antennas are small (40 cm 
long) and are not likely to be damaged while the buoy is hand¬ 
led. Both the liqht and the radio have their own batteries; 

both can operate for two days. 

In typical operation, the buoy is started, and placed in 
the water. Some time later it starts and data is recorded. It 
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is difficult to stand around hoping all is well. To aleviate 
cistrcss in the operators, the buoy has a small (100 m watt) 
citizen-baid iadio transmitter that .ransmitts everything it 
hears. The operators, several hundred yards away cdi hear the 
gyros turn on and spin up, and can hear the tape recorder step 
each time data is recorded. This verifies that all is well, 
and the operators are relieved. 

The general layout of the buoy hull is shown in the cut¬ 
away view (frontispiece) and figure 4. The control logic, ver¬ 
tical gyro, inverter, and transmitter are mounted together on a 
frame (figure 5). The data recorder (figure 2) mounts directly 
on top of this. The entire assembly is housed in a watertight 
aluminum can. The light and transmitting antenna are mounted 
on top of this. The north-seeking gyro (in a waterproof poly¬ 
vinylchloride container) and the batteries mount on the buoy 
hull. Underwater type connectors connect them to the lid of 
the can. The container is 52 cm in diameter, 40 cm high, and 
weighs 60 kgm. 

The hull is made of two, 1.5 meter (5 feet) diameter 
sheets of 1" marine plywood. Ten centimeters (4 in) of closed¬ 
cell styrofoam is sandwiched between them for floatation, and 
the three bolt together to form the hull. The aluminum can and 
gyro bolt into holes in the top piece of wood. The batteries 
are held below three hatch covers. All wiring is routed between 
the plywood sheets and protected. A small wind vane provides 
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further protection for the compass gyro. The edge of the hull 
has a soft neoprene bumper glued to the styrofoam. 

The entire buov assembly bolts together. The aluminum 
can is easily remo' ed and carried into a protected laboratory 
for repairs, calil ration, and setting switches for turning on 
the electronics. The electronic rack may be unplugged and 
removed. All electronic components, gyrco, etc., plug ir and 
can be removed for repair. Connections to the data recorder 
are by screws or. a terminal strip. All internal wiring is 
laced together. Thus the entire assembly may be quickly and 
easily disassembled for repairs. In normal operation, the unit 
can be switched on, the lid to the aluminum can bolted on, the 
can bolted into the hull, and the buoy launched, all in 20 
minutes, by three people. 

A three-point bridle of spliced 3/4" nylon rope is attached 
to the buoy for launching. Usually, 30 meters of 3/4" polypro- 
pvlene rope is attached to the bridle. Vihen the buoy is in 
the water, this line floats on the surface upwind of the buoy. 
The heavy ]ine is necessary to handle the buoy in rough seas. 

A smaller line (1/2") parted during one recovery. 

The aluminum can ana wind vane are painted bright yellow 
(epoxy coating). The hull is bright red. English and Japanese 
lettering on the vane identifies the buoy. External fittings 
are stainless steel, external wiring is neoprene jacketed. This 
greatly reduces corrosion and electrical shorts to sea water. 


5. DATA ANALYSIS 


The buoy data are analysed on a Burroughs B6700 computer. 

A complete listing of the program is produced in the appendix. 
Although it appe rs long and complicated, it is, in fact, 
straightforward. Processing of data is quick, typically one 
hour of buoy data is analysed in 11 minutes at a cost of about 
$25. The following paragraphs outline the procedure, shown 
graphically in figure 6. 

Data from the four sensors {acceleration, two tilts, head¬ 
ing) are written on magnetic tape as six-bit characters. The 
first 24 bits of each record is a time word. The procedure 
TAPESTART positions the tape at the record to be read, and notes 
any read errors in records passed over. The procedure READTAPE 
transfers into the computer the bits written in a particular 
record, converts them into time series with appropriate units, 
removes the mean, and notes any apparent errors. Procedure 
PRDATE decodes the first 24 bits into time and date. The four 
time series are trai ^formed into the frequency domain (using 
the Fast Fourier Transform) by the computer's library routines 
FFT and BITRV2. 

The frequency data are multiplied by the appropriate 
inverse filter function L -1 calculated from (4.1). Some data 
were recorded with sample-and-hold modules installed in the data 
recorder. In this case, all channels were sampled simultaneous¬ 
ly. Other data were recorded without the modules and the chan- 































nels were sampled sequentially. This is equivalent to applying 
the linear filter: 

L = cos nut + jsin nut (5.1) 

where n = 0,1,2,3 is the order of sampling, and t is the time 
delay between samples.. When appropriate, the inverse of this 
filter was applied. Finally, the tilt and heading were trans¬ 
formed back into the time domain. 

At this point, the slopes are relative to the buoy, and 
must be transformed into a geographic coordinate system. Let 
r be the rotation about the outer gimbal axis, p the rotation 
about the inner gimbal axis, and °T the angle measured clock¬ 
wise from North to the inner gimbal axis in the coordinate 
system shown in figure la. The coordinate transformation is 

; 2 = (-sin°T/cos p)tan r + cos°T tan p 

= (-cos°T/cos p)tan r - sin°T tan p (5.2) 

These equations follow from Saenger's (1969) equations (9.6) 
for ip = 27C° and ip = 0°, and noting m = -°T + 90°. 

The slope series are transformed back into the fiequency 
domain and the appropriate Co- and Quadrature-spectra are com¬ 
puted from the Fourier coefficients of the data using: 

2C ij - a i a j - b i b j 

2Qij ~ a i^j "*■ d jt>i 


(5.3) 


where a^, are the cosine and sine transformation of the 
first series, a^, are those of the second series. Note 

that the sign of is arbitrary. Our definition agrees 
with (2.9) , but is .n disagreement with the convention used 
by some workers. 

The spectral quantities (2.9) are used to calculate 
various parameters of the ocean-wave directional spectrum. 

The mean wave direction (at each frequency) is: 

tan fc> 0 = Nqi/Niq (5.4) 

The root-mean-square beamwidth can be calculated several ways; 
we use 

tan 2y * (rrl ) V2 < 5 - 5 > 

where 

r2 = t (N 20- N 02 )2 + ^i^/Noo (5.6) 

This agrees with equation (22) in Longuet-Higgins, Cartwright 
and Smith (1963) . 

6. BUOY RESPONSE TO WAVES 

We have estimated the response of the buoy to waves by 
observing the internal consistency of data recorded at sea and 
by calibration in a wave tank. In general, the buoy responds 
as a damped harmonic oscillator with a 1/2 power point near 


0.7 Hz. The accuracy of the data meets the design criteria 
specified in section 3. 

Before describing the calibration of the buoy, we first 
give a gross overview of some typical data. They were obtained 
off Monterey, California in a 10 m/sec wind when the signifi¬ 
cant wave height (4<^ 2 > 1 / 2 ) was 2.5m. A short section of the 
digitized data from the fo ir transducers (figure 7) indicate 
the general range of these variables. The root-mean-square 
acceleration was 10.3% g, and the slope was 6.8°. This is in 
very good agreement with the values of 10.7% g and 6.1° predic¬ 
ted by eqs. (3.4) and (3.5) using Q = 0.7 and v = 10 n;/sec. 

The first time the buoy was operated at sea we sampled the 
transducers at the maximum possible rate to determine the opti¬ 
mum rate for later operations. The Nyquist frequency (maximum 
frequency of the digital spectrum) was 6.25 Hz. The accelera¬ 
tion spectrum recorded on this occasion is shown in figure 8. 

The noise at low frequencies has a spectral power of around 
10 cm 2 sec -4 Hz“ 1 . The high-frequency noise is considerably lower, 
and the Nyquist frequency could be reduced considerably. Conse¬ 
quently, all later data were recorded with a Nyquist frequency 
at 1.5625 Hz (each signal was sampled 3.125 times per second). 

The spectra derived from the Monterey data are plotted in 
figure 9. The Nyquist frequency is 1.5625 Hz, the resolution 
is 0.0015 Hz. The figure gives a good indication of the signal/ 
noise ratio of the data. The acceleration values are approxi- 
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niately 1000 times the low-frequency noise, and the low and high 
frequency ends of the spectrum are identical to those in figure 
8. The slope spectra are approximately 100 times larger than 
the low-frequency noise. The slight rise in energy near zero 
frequency is due to drift in the sample-and-hold modules, and 
is the reason for not using them on subsequent days. The buoy 
heading swings slowly with time, so its spectrum is large at 
low frequencies. The small spikes in the spectrum are due to 
aliased 400 Hz noise. 

The acceleration noise values require some comment. The 
noise is 30 times that expected from least count noise (which 
itself is about equal to the noise from the transducer). It is 
also about 30 times larger than the noise measured with the 
instrument on land (no motion). It is probably due to the 
slightly non-linear way the buoy responds to the motion of the 
sea, and is not significant for our work. For example, a sinu¬ 
soidal wave with an amplitude of 0.8 cm at the lowest frequency 
of interest (0.06 Hz) is detectable with a signal/noise ratio 
of one. The contribution to this same frequency band by a 
fully developed sea is 100 times larger. 

The buoy was calibrated by observing its response to 
sinusoidal waves in a tank 2.6 meters wide, 2 meters deep, and 
30 meters long. This facility is just barely adequate; fortu¬ 
nately, the best data was obtained near the buoy cut-off fre¬ 
quency. Lower frequency waves were too long for the shallow 
tank, and shorter waves were not two dimensional. The observed 
response is plotted in figure 10, together with the 
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response of a 1/8 ellipsoidal hull calculated by Kim (1966) . 
Although this hull is not very similar to the buoy in plan view, 
it has nearly the same cross-section along one axis, and it is 
of shallow draft. The wavenumber k of the incident waves was 
calculated from the wave frequency w using w 2 = gk tank kh, 
where h is the depth of the water. The figure indicates the 
model accurately predicts the 1/2 power point in the buoy's 
response, and the ratio of pitch to heave response at low fre¬ 
quencies. This ratio diverges from the observed values at high 
frequencies. 

We estimate the internal consistency of the data by 
observing how well eqs. (2.10) and (2.11) are met. In general, 
C 12' C 13' and ^23 are zero with i n the statistical fluctuations 
of the data. Six hours of data recorded in a trade-wind sea, 
when averaged to give spectral estimates with over 2000 degrees 
of freedom, yields values for these quantities that are less 
than 5% of the leading terms (C^) . To evaluate (2.11) wp have 
averaged together 17.23 hours of data and plotted it in figure 
11. The peak near 0.9 Hz is predicted by Kim (cf. figure 10), 
but the slight rise above unity (near 0.5 Hz) is unexplained. 
Perhaps it is due to the additional moment of inertial of the 
high instrument can. 

7. CONCLUSION 

The wave ouoy has worked well. Typically we work near 
shore. The buoy is taken 20-30 km offshore, data is collected 
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for three hours, and we return, all in one day. In one year, 
it has operated at sea on nine different days in four different 
parts of the world and recorded more than 30 hrs of data. Dur¬ 
ing this time there was only one major failure. On the last 
day of operation the unit stopped because of low battery volt¬ 
age. Inspection revealed sea water had corroded away most of 
the terminals ana discharaed the batteries. 

The data recorder has also worked well. We usually record 
data in 10 minute segments (records) of 2048 scans. Out of 
over 200 records recorded, 5% (11 records) have errors. All 
but two are due to the tape recorder adding an extra 6 bit char¬ 
acter in the middle of a record. This is not serious because 
the data can be recovered with a little extra computer work. 

We prefer to ignore these records. Once the clock failed to 
act properly. No data was lost but the time word was incorrect. 
Once the record was' too long. We suspect these errors were due 
to noise introduced into the logic circuits by strong radar 
signals from ships used to deploy the buoy. Certainly, the 
errors were less frequent when the ship did not have a radar or 
was not close to the buoy. 

The buoy has withstood rough handling. It has been shipped 
from San Diego to Monterey and back, to San Clemente Island and 
back, to hawaii and Wake Island. The waterproof connectors have 
been damaged by inexperienced crane operators on the ships, and 
the radio antennas have been bent. These items are easily 


repaired. The internal components, which are harder to repair, 
have not been damaged. 

The buoy has been deployed by helicopter ana from a variety 
of very small ships with little difficulty. Two men holding 
lines attached to the bridle mounting points can steady its 
motion on a rolling ship. It can be lifted by small hydraulic 
hoists or cranes found on many small workboats, and an inexpen¬ 
sive two-man helicopter can carry it to sea. 

Our experience with the buoy shows ways it can be improved. 
The instrument can could be considerably smaller. It was orig¬ 
inally designed to hold the compass, but iron in the inverter 
caased interference. If only three hours of data are recorded 
at one time, the much cheaper low-power, cassette type data 
recorders could be used. This would make changing tapes easier, 
and would further reduce the size of the instrument can. Thus the 
buoy \.ould be more disc shaped, and should have a better response 
to waves. 

In conclusion, we can say the buoy has met the design cri¬ 
teria specified before it was built, and has provided the ocean 
wave spectra we need in our research program. 
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APPENDIX 

This section contains many of the details necessary to 
operate the buoy and to analyse the data from it. We have 
included detailed circuit schematics (figs. A1, A2, A3), pro¬ 
gram listings, and a list of major components and their source 
of supply. 
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Listing of Computer Program used to Reduce Buoy Data 


i LC1> 

V 

i program, wavebgoy provides a means of computing estimates of 

» Thk DIRECTIONAL SPECTRUM OF OCEAN «AVE5 FROM DATA TAPES 

Gr. NEnATEU MY 1IE S. 1. U. SURF ACl-F U LL'J W I N 0 WAVEBJJY. 

% 

THt CARO INPJT PARAMETERS TO BE SUPPLIED ARE: CARL 1, THE TAPE 
•, RtEL NAME (0 CHARACTERS); CARD 2, THE T IfE-SER IES LENCTH, 2«*MEX. 

t FOLLOWED BY THF NY'.’UIST FREQUENCY IN HEPT2. CARO 3 CONTAINS FIVE 

. INTFCEP CONSTANTS VIZ..THE FILE NO. JF THE DATA, I HE FIRST RECORD 

OF THE DATA TO BE USED, THE NUMBER OF RECORDS TO BE AVERAGED, 

, AN I NY EGER VARABLL,“CHANNELS",WHICF SELECTS WHICH OF THE FOUR 

\ DATA CHANNELS ARE TO BF USED. IF ALL CHANNELS ARE REQUIRED,AS IN 
E COMMUTING OIKECTIONAL SP-CTRA,THEH SET CHANNELS:*1111 (BINARY) * 

A 15 (DECIMAL). FOR INSTANCE.IF ONLY CHANNEL 1 POWER SFECTRUM (VERT, 

t ACCELERATION CINVENTED TO SURFACE lLEVATION) THEN SET =1000 
\ (dINARY) * 8 (DECIMAL INPUT VALJE1.“AVEFJ"IS THE NJMBLk OF KREJ. 

A BANLS AVERfGEl OVER AFTER THE AVERAGING OVER SPECTRA I'OM 

t INDIVIDUAL RECORDS HAS d t E ,N DONE. THE NEAT CARD READ l.iPJTS: 

\ A BOOIEAN VARIABLE VAV" «HICH.IE SET TRUE,DOES A LOCFITHMIC TYPE 

i AVERAGING OVER SPECTRA^ BANDS (CAUTION:IF LAV*TRUE,SET AYEEU*!). 

* A BOOLEAN VARIABLE “S Ahl N", WHICH IE SET TRUE MEANS THAT T,. r 

t s i M J L •' AMD HOLD MODULLS OF THt WAVEdJOY WERE OPERATIONAL AT THE 

TIME OF LATA COLLECTION; THE MAGNETIC COMPASS DEVIATION FkUf TRUE 
} NORTH ANO EITHER “E" OR "W" CHARACTER FOLLOWING; "TLIM" IS 

i Tht OKEATEST WAVE PtRIOD 10 bE INCLODEI IN FINDING THE VARIANCE;IE 

bLANK THIS IS AUTOMATICALLY SET TO 2* * ( M E X-1 ) / NYU J I S T, Trlt LONGEST 
> PtRlOL FOR 1 HE TIME SERIES LENGTH. UNALLY, THERE IS A CARD 

A WHICH CONTAINS SIX HCOLEAN VARIABLES SFECIEYING WHICH TAbLES ARE 

C TO dE PRINTED OUT AND WHICH TYPt(S) OF SPECTRA ARE TO BE PJNCHED. 

» 

; THIS PROGRAM UTILIZES Tilt %StT IN3ALLA T ION" CONTROLL CARD 

4 TN ORDER TO ACCESS THE SYSTEM INTPINSICS: FFT,EITRV l ,F)TR,SINCOS. 

* JOE JOY , TAOl HANSEN 

v BOB STE»ART 

APRIL, 1873. 

*, » » » 1 I M t I ! » Y I M I M I I I 1 M M 1 I I » I M I 


nno : finnn : r 

BLOCK IS 3ECMENT OflODt 
nr 2 : 0 nnti : r 
002:0000 ; r 

002 : POCO;0 
002:0000:0 
002:0010:11 
002 : ( ( iOf';fi 
002:0000:0 
002 : 0000 : l ‘ 
002 : 0000:0 
002 : 0010:1 
002:0000:0 
102 : 0010:1 
I 0 2 : r 0 0 ft : 11 
00 2:00 ) 0:1 
00250000:0 
002 : 0010;0 
002 : 0010:1 
002: 00(>0.: 0 
1102 : 0010:1 
002 : 0000:0 
002 : 0011:1 
002 : OP00 ; i 
012 : 0110:1 
002 : 0000:0 
002 : 0010 : i 
002 : 0000:0 
002 : 0710 : ) 
002 : 0000:0 
no 2 :norn : o 
002 : 0011:0 
002 : 0000:0 
112 : 0110:1 
002 : 0000:0 
012:0010:0 
002 : 0000:0 
002 : 0010 : i 
002 : 00 ( 0:0 
002 : 0010:1 
002 : 0000:0 


HLE PRIN TEk (KIND * 6, HAXRECSIZE = 22); 

FILE CARD (KIND = 9, PAXRECSIZE = 14); 

PILE PUNChEk (KIND * PUNCH, flAXRECSIZE * 14); 


LOOLEAN EKK; 

IMIECPR CHANNELS. FIRSTREC. J, K, LASTREC, HEX. N. NHl, NCHANNEL, NR, 
NULE, NEC, NWORDS, NREC, TAPENAME, AVEFQ; 

PEAL FSTEP, OMEC, SCALE, VARN, X, XNYQUIST; 

REAu ARRAY CAL, HEAN(1:4); 

I EH Pt PI * 3.1415926535897932W; 

A 

A NOTE: 3777 OCTAL * 2047 DECIMAL * *5.0 VOLTS. 4000 CCTAL*2048 

FILL CAL[* ) WITH . 4 78515625 A ( 983 CN/SEC ** 2)/2046. 
,.001533980768 A PI/2046. 

,.000383495197 A (PI/4)/2046. 

,. 000 38 3495197 ;A (PI/4 )/2048 . 

PEA I (CARD, <.'b> .TAPES AHE); 

READ(CAU),< 14, F10.4>, HEX, XNYfaUIST); 

KEADICARD, < 51 3 >, NFILE, FIRSTREC, NkEC, CHANNELS. AVEP fa); 

1) AVSFJ LEJ 0. THEN AVEFQ:*1; 

V>Rni(fkINTER,<" TAPE REEL*", A6>, T APENAHF.); 

WklTK PkINTEk.<“PEX*",l3,X5,"NY0UIST*".P8.4,“ HERTZ.">,HEX.XRYQUIST); 
NWORDS :* 0 s. l[HEX:3:l); A 2 ** HEX 
N :* RWOPDS DIV 2 ♦ 1; 

RRl :* N -1; SCALE:*2./NW0RDS**2; 

A 

WRITE! PRINTER.<"PILE N0.*“,I3," BEGINNING WITH RECORD NO. ",I3, 

" TOTAL NO. OF kECORDS READ*",13, 

" CHANNELS*", I 3>,NFILE,FIRSTREC,NR EC, CHANNELS); 

A 

A 

A 

A 

BEGIN A <<<<<<<<<<<<<A HAIN PROCRAH BLOCK BEGINS HERE. 

A 

A 

A 

B3DLEAN SAHIN, TABLE1, TABLE2, TABLE3, PUNCHX3, PUNCHPS. PJNCHOS, LAV; 

1NTECEk EAST.E OkW , ENDPT,FkEEDOM,VALAVEFQ, WEST; 

REAL C, Ci1. C33. C44, C34, CCSAZ, COSPITCH. K2, Kl, HAGDEV, fal3, J14, 
kk, SINAZ, TLIH, TANPITCH, TANROLL, Y. ZETA, Z; 

PEAL NZERO, DER0H1, DEN0H2, JACOBK, Kll, K22, N01, N10, N20, N02, Nil, 
N12, N21. N22, C14, C13, Q34; 

REAL ARRAY 4. I [1:4, 0:NW0RDS), PSA, A, B [1:4, 0 ; .NH1), 

CSN. CSW. CNW, CO, FREQ, 0HEC4, PER, PS, faSN, CSW, faN*, 

I I[0:NH1 ), AVEPlR, A3, Al, A2, Bl, B2, BETA, CREST[n ; 1 t 
ENTIER(NP1/AVEFQ)); 

PEAL ARRAY NAS, NA7, NAB, RLD.NAl, NA2, NA3, NA4, NA5 , ANG1 , A.NC2, RNEW, N.NN, 

XALPH [ 3 : 1 ENTIERI NH1/AVEFfa ) ), JACOB,K0NE,KTW0[0:NH1J; 
UPECT FILE BUOYTAPE (KIND * 13, LKBELTYFE = OHITTEDEOF, DENSITY * LOW, 

PARITY * 0, HAXREC3IZE = NWORDS ♦ 1); 


I'lHECT ARRAY TAPEBJPI0:NW0kDS|; 

UDLL ADD Ik, AGAIN,FXIT.FINI,I.OOPER,.NEXT,PUN1,PUN2,PUN3,TAb2,TAB3; 

% 

<<<<<<< <<<<<<<<<<<< 

P3«r.Ar p 5( " N ",X5,"PERIOD",Xd."A(U",X12. " Al * ", X1 2," A 2 *" , X12," Al *", 

IKL)’.",X5,"PfiS i/i.GPEES",/ ); 

1AT HI It. <2,r 3. 3, X l ,5(U2.4, X3 ), 2< K7.W , X7 ) ) ; 
rOPHT EMC* ANG.", XI. "PERIOD", X3. "ACCEL.", X4. "SURF ACE", X3,"N-SuOFE", 

X3,"t-SLOPE", an.X t,"QJAD-S:C , X4 ,"C()-SE", X4, "Q0AD-3E", 
X4,"CO-NE",X4,"faUAI-Nl".X2,"WAVES R F A H" ); 

) 3RNAT Fl'll Kl H,"FK3H", X5,"WDT4" ); 

FOPHPT UK" FkF.u.".X5,"SEC.",X3,"'H2/SEC4",X3, M (H**2",X4,2("PER hEkTZ 1 
).4(" CN./HZ. "),2( "PER HERTZ ")," JEG. (T)"); 

I0P7AT H2( E6.2,F6.2,10E10.2,21X2,F6.1)); 

I3R1AT Fill" AVEHACED DIRECTIONAL COEFFICIENTS FOk .969 SEC. < T < ". 

F 7.1," SEC.",//); 

flklAT F14IX32,"WAVES FR01".X5."BEAT WIDTH"); 
rUPMAT ) 1 5( X18,"PEP HEkTZ",X2,"PER HZ.",/); 

FOkTAT )UI "ANC. )kE5.",x5,"»ni/Non",x5,"Nin/NOn“,x5,"N20/NriO".x5. 

"Nn2/Nn",X5,"Nll/.Nno M ,x5,"ALPHAl",X3,"ALPr(A2',X3,"RNEW", 
Xt,"ROLL",X5,"DIXP",X3,"ALPHA",/); 

IJRIAr Fl/(El0.4.x9,S(ES.2,Xi),2(F6.1,X2).2(Ei.2,Xl).f 0 ./,X2.t9.2); 
.•UVA. Hbt Arc. KF w .",x5,";ii2/Nor",x5,"N2l/Nnn",x5,"N22/NTr ", <t, 

"NCHECK", /); * * 

FJiKAT Ft MElfl.4,XN.4(29.:,xt'); 

< << <^<<<<<<<<<<<<<< 


002 : 0000:0 
DATA IS 0005 l3nC 
002:7000:0 
DATA IS 0003 LONG 
002 : 0000:0 
DATA IS 0003 LONG 
002 : 0000:0 
002:0000:1 
002:0000:0 
002 : 0000:0 
002:0030:0 
002:0002:1 
002:0002:1 
002:0002:1 
002:0002:1 
002:0001:5 
002:0003:5 
012:0003:5 
002:0006:3 
002:0014:5 
002:00*5:2 
002:0038:5 
002:DC 3D:4 
002:07 4B: , 
002:03 36:1 
002:0051 :1 
00 2:005F:0 
002:0063:2 
002:03,3:2 
002:0065:1 
002:0035:1 
002:0077:1 
002:00/7:1 
002:0077:1 
002:0077:1 
002:0077:1 
002:00/7:1 
002:0077:1 
002:00/7:1 
002:0077:1 
002:0077:1 

BLOCK IS SEGNiNT 00015 
005:0000:0 
005:0000:3 
005:0000:0 
005:0000:1 
005:0000:0 
005:0000:1 
005:7008:3 

005:0018:3 
005:0011:0 
005:0319:1 
005:0119 :1 
005:01 *9:1 
035:0329:4 
DATA IS OP25 LONG 
DATA IS 00:13 LONG 
OO5:01*C:4 
005:0030:4 
005:1330:4 
005:7030:4 
005:0330:1 
005:07 jO;4 
005:0330:4 
005:7r30:4 
018:0330:, 
005:07 30;4 
005:03jO;4 
005:0030:, 
005:0331:4 
005:0(1 jO: 4 
003 : 0130:4 
005:0730:4 
095:0330:4 
005:0730:4 
005:0030:4 
005:03 37>: 4 
005:0730:4 
005:3330:4 
005:0730 :4 
005:0730:4 
005 : 0130:1 
005:0030-4 


A. 




I 


K Rtc PLUrL 

t 

TKKEST'AR i 

(lIRECTPiLt, OB, 

RS 

. SK, 

SR, lRR >; 

V AIDl no. 

si, sr; 

INTEGER RS, 

SK 

, SR; 

BOOLEAN lRR; 

Lii-ECT KILL IIRECTF ILL; 


DIRECT 

ARRAY Lb 10); 

>.4444 

t 4 4 4 4 

t 4 4 4 t 4 4 t 

4 4 

4 4 4 

4 4 t 4 4 t t t 4 4 


t 

» ,'HIj PROCEDURE IS FOR k-.ADING A HELD DATE LOGGER TAPE BJf 1EANS 

* OF DIKKCT I/O. THE UK IS ADVANCFD TO A SPECIFIED KILL GIVEN 

» hi "SF" AN L TI.E1I TO A PECcRD GIVEN BY “SR".THE EXPECTED NUKEER 

\ OK »G41j EACi RECORD CONTAINS IS "RS". IHt ACTJAL NUflbEK OF WORDS 

A IkANSFEREL KKOn EAcll RECORD AM Thl READ ERROR TYPt AkE PRINTED 

JOT FOR EACH KILE AS WELL AS ME SUIBER OF RECORDS PER KILE, 
i PRJCEO )RE TAPKGTArT REQUIRES KJREWAkD DECLARATION CF DIRECT FILE 

I AM DIRECT ARRAY DB. 

*4*44444444444444444444444444444444 

■> C. J I N 

INTEGER EOKCT. RC; 

KILL PRINT IKINL * PRINTER, f.AXKECSUE * 22); 


on5 : nnjn:<; 
n.ia.-mjn:! 
nob:OC 30;A 
0Ob:0 T 30; 4 
005:0130:4 
Of'b: DC 3D: 4 
005:01 j‘I: . 
OC5:0C30:4 
00b:0r 30:1 
005:0130:4 
005:0fjO : 4 

015:0131:1 
005:0131:4 
00 5:Of 30:4 
003:0131:. 
003:01jO:4 
103:0131;; 

3 (05:0(30:4 

TAPESTART IS SEGMENT 00107 

0 (i 7 :f onn : 0 

DATA IS 005A LUNG 
DATA IS 000 3 LONG 


KORffT SKRC ("KILE”, 13. " RICORI", 13, " 
"EkRJR NO.". 12); 

WKITLIpRINT(SPACE 2),<”SPACE TO PILE ", 
tPR := KPLSE; 4 NO ERROR INDICATION. 

«C :* 1; 

LLI LV := 1; 

EOKC: LJS SK DO 

ELGIN 

READ I JIRtCTHLK, RS, L'd[* |); 

A a IT (l B [ * 1) ; 

IK oa[*j.11L3K THEN 
btGIN 

CiiOSEl DIA EcTK ILt,'); 

WRITE ( F RINC(SPACE 21. <“FOK ENCI 
13, " rECJRO(S)">, EjPCT, RC); 
ECFCT :**♦!; 

DIRESTPIL_.1FEN :* TRUE; 
r C : J 0; 


CONTAINS". 15. 
12. " RECORD ' 


END ELSE 
BEGIN 
C 


» 1; 


t NO 

end; 

kC :- 1 i 

WHILE RC LSS SR L0 
BECIN 

Rt.V (DIRECT FILL, RS, U B (' )); 
W AIT ( D B [ * ] ); 

aKU'EI PRINT,S'KAC, t'JFCT, RC, I 
kC :« * ♦ I; 
tvu; 

P..I TAPESTART; 


WORDS. ", 


007:0000:0 
007 : 3110 : 0 

12>, SF.Sk); 


007:0000:0 

007:0010:5 

117:0111:4 

007:001^:2 

017:0013:0 

007:OOiJ:b 


4 

117:0114:2 
007 : Of la: ^ 
0 fi 7:001 a; 1 
007:DO 1C:0 


5 

117:011C:3 

", 12, " OF", 


0O7:001t:0 
007:0121:3 
007:f'P2t:5 
.107 :0030:1 
007 :(.031:5 
007:0032:3 


5 

007 : 00 j2 : 3 


5 

007:0033:1 

ioerkoriype ); 


007:0034:2 
10 7 :00 4 b:1 


5 

007:or 4c: 1 


4 

007:0148:4 
0 0 7: T C 4 C: 2 
007:0040:1 


4 

00 7 :01 4D:» 
007:0051:4 

RRuR TYPE); 


017:0033:2 
007 :OObA: 1 
007:0058:3 


4 

007:00bC;O 


TAPESTARTt007) 
3 


IS 007d LIN 
0 05:0 f JO:4 




IkOCFDURE PKLATL ( DAI ); 

f 

t 

VhI Lt DaT; IhUGER DAT; 

V 

it M l M M M » M » » » M » M U U U > > M M l 

t 

* T ‘IIS PPOCLLi'W; RECODI S TH( DAT' AND TI"t INF ONPATI ON CONTAINED 
\ It; THE LLMItlC 2* PITS Cl FfCH TAPF prCOPD AND PRINTS THIS OUT. 

V i t t t t < i < t t t < t l V t t t « t t t t t < i t t t t < t t t 

i) t J I N 

INTEGER ARRAY rOMHS[l:I2|; 

VALUt ARRAY MONTH I "JAN", "r EH", " 1 AH", "APR", "MAY", "JJN", "JUL", 

"AJG". "SEF“, "OCT", "NOV". "DEC"); 

HTcCEr YEAR, DAYS. HOUR. PI IN 1^, MINI, NH; 

UELL INI. EXIT; 

FILL ’iONHSl*) WITH 31. 20. 31. 30. 31. 30. 31. 31, 30. 31, 30, 31; 
CJTHE'.r JAT (DATE AND TINE) IS IN THE FORK: 

Y ID lLCD ULD IH HHHh YrtP. MMPM Y-YEAR, D-DAY, H-HOUR, r - PI I (4 
IN 9INARY COJtD JnCiiiAL DIGITS; 

YEAR :* IF CAT.[23:1| EwL 1 ThLN 7j ELSE 72 ; 

*AYS :* DAT. [2<: 2) * ion * tAT.[2'i : 4| * 10 ♦ DAT.(ifc:4 ); 

I,DDK :* DAT.[12:2] * 1 r ♦ IAT.[10:4|; 

MNlc :* tAT.! t>: 3 ] ; 

TINl : = DAT. 13:4]; 

IF YEAR POD 4 ECL 0 THEN P ONTHS[2 ) :« 29; 

F3R Nn :* 1 STEP 1 UNTIL 12 DO 

IF CAYS LE V MONTHS ( NP.) THEN GO FI.D ELSE LAYS PAYS -MONTHS [ N PI)» 

» KITE(PHINTtF,<"ERRUR IN PRDATE: DAYS TOO LARGE"'); 

CO tXIT; 

END: WRITE! PRINTER. <"DATE: ", 12, XI, h3, Xl, 12, "; TIME: ", 12, 
II, 11 >, DAYS. FIUliTH [ NP-1J, YEAR. HOUR, niNlO, MINI); 
nimtlPRlNTEK.c/*) ; 

EXIT: FND FRCATt; 


005:0331:4 
onb:Pf jO: *, 
nnb:PC30:4 
005:0330:4 
005:Of 30:4 
005:0130:, 
005:0030: 4 
"05 : 0 13 3 :4 
00b: (-r 3n : ^ 
nn5:0r30:4 
015:033 ':4 
005:0030:4 
005:0330:4 
0.05:0T30;4 

PR LATE IS SEGMENT 00005 
039:0030:0 
DATA IS 001L LONG 
0n9:irnn;0 
DATA IS OOOC LONG 
009:0000 : 3 
009:0000:0 
009:0000:0 
009:0014:2 
009:0004:2 
10 3:OO14 : £ 
009:POOH :2 
004:0017:5 
005:PO OC:2 
005 : 0001 : 2 
004:0011:5 
009:TO 12: t 
004:0015:1 
009:0015:5 
OO9:0OlC:E 
009:0021 :5 
n09:0022:2 
009:0024:1 
009:0039:1 
005:0040:4 

PRDAIE(009 ) IS 0043 L0N3 
005:0030:4 


■,<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 


3 


PRIZEUUKt KEADTAPE! NFILE. RtC. CHANNELS); 
>. 
i 

VALUE NFILE. NEC. CHANNELS; 

INTEGER NFILE. hEZ. CHANNELS; 


% 4 \ 

BECIN 

I 


HmnmnntumnnnmtiMn 

THIS rROCEDURF READS ONE RECORD/CALL BY MEANS OF DIRECT I/O. 

THE HIGH ORDER 24 BITS OF THE FIRST WORD. WHICH CONTAIN THE DATE 
AND TIME INFORMATION.ARE ISOLATED INTO THE VARIABLE ‘DATE*. 

ThE FOD.' CHANNELS OF DATA ARE SEPARATED AND CONVERTED INTO THEIR 
APPROPRIATE UMTS. THE P.EAN VALUE OF EACH TIME-SERIES IS FOUND 
AND SJ3TRACTEO OFF. VARIOUS ERROR DIAGNOSTICS ARE PRINTED OUT. 
RtADTAPE REQUIRES FOREWARD DECLARATION OF TAPEFILE A DIRECT ARRAY. 


t % 




NTEGEW date, errct. ioerr. j. l. limit; 


REAL CF. XMEAN; 

BOOLEAN ARRAY S WCH (1:4); 

LABEL CCNT, ERRLP., LEAVE; 
OEFINt SETFIELD! SUBC. L3F1 * 
LNTIL J :* * * I CEQ NWORDSA; 


OU R(L.J) TAPEBJF|SiiaC). (L8F:12) 


DIRECT 


i/o; 


" RECORD". 13. 
REC. TAPEBUF(* 1 


' CONTAINS". IS, 
IOWORDS. IOERR) 


IF 


ERR:* FALSE; 

Z01MENT READ THE NEXT RECORD USING 
READ!BUOYTAPE. NWORDS*l. TAPEBUF(*] 

WAIT! TAP EB JF [ * ] 1 ; 

IOERR := TAPEL*UF[ * )- IOEhKCRTYPt; 

WRITE!FRINIFA[SPALE ?]); 

WRITE! PRINTER, <"FILE NO.". 13. 

" WORDS ERROR NO.". I2>, NFILE. 

CASE IGEhK OF BECIN 
; t 0 NO EHaJR; 

BECIN WHITE!PRINTER,<"NOT READY">); CO TO ERROR ; cND; 

TAPEBUFI* 1.IOWOKLS NEQ NWORDS ♦ 1 THEN GO ERROR ELSE WRITE!PRINTER, 
<"PANITY ERROR. RECORD WILL BE PROCESS ED.">1; » 2 PARITY. 

BECIN WHITE! PRINTER,<"REWIHuING">); CO TO ERROR; END; 

WRITE(PRINTER.<"TESC. ERROR">); GO TO ERROR; END; 

WRITE! PRINTER,<"END/BEC. TAPE">),* CO EXIT; END; 
WRITE(PRINTER,<"FND OF HLE">); CO EXIT; END; 

•RITE!PRINTER.<"ERROR A7">); CO EXIT; END; 

WRITE! PRINTER,<"ERROR * 8"'-l; CO ERROR; END ; 

WRITE! PRINTER.<"SHORT RECDRD ".FlO.3," CHARACTERS">, 


BEGIN 

BEGIN 

BEGIN 

BEGIN 

BECIN 

BEGIN 


TAPEBUF!*).IOCHARACTCKS); CO TO ERROR; 


END; 


END; 

l : * 1; 

DO 3 'WZ H [ L ) :* BOOLEAN! CHANNELS. ( 4-L: l 1) UNTIL L :« 
COMMENT DJMP HEAOER INFORMATION ON THE LINE PRINTER; 
CONTs 

DATE :* 0 fc TAPEBUF!0 ] [ 23: 47 :24 ); 

P RCATE(DATE); 

FOR L :* I STEP I UNTIL 4 DO 
IF SWCH[L] THEN btCIN 


♦ 1 CTR 4; 


3 

READTAPE IS 


4 

5 


b 

b 

b 

b 

b 

b 

b 

b 

4 


PO b:0C30:4 
70b:033P:l 
005:0030:4 
OPb:03 3P: 4 
OPb:0(30:4 
00b:03J 0:1 
005:01.40:4 
035:0130:4 
005:0330:4 
0Cb:0C30: a 
00b:0030 :a 
00b:07 30:4 
005:0030:4 
00 5:CC 30:4 
005:0130:4 
00 b: 0 f30:4 
005:0030:4 
00^:0731:4 
005:0030:4 
005 : 0130:1 
SEGMENT 0010 
OOC:0000:0 
IOC:0030:1 
OOC:0000;0 
onc:OOin : .i 
OOC:0000:0 
00C:00in;0 
OOC:0000:0 
n0C:0130:4 
OOC:0000:4 
03C :003 4:b 
000:0006:3 
nrc:nong : i 
OOC : 003 E:3 
00C:00in : 2 
OOC:00/4:1 
000:0026:3 
OOC:032b:3 

OOC:002D:0 

OOC :0031:3 
0OC:0036:2 
OOC:00 3C:b 
000:0043:2 
OOC:004b:0 
OOC:0052:4 
OOC: 005A:2 
OOC:Onb3:b 
OOC: 0063: 3 
OOC:007 1 :1 
OOC:007b: 3 
OOC: 007 7 :1 
OOC:007C:3 
OOC: 007C: 3 
OOC: 007 C:3 
OOC:C07t:3 
OOC: 03 7 F:i 
OOC:OObO:/ 


J := f; 

CASE L-l OK BEGIN 


GETFIELLIJ.23 >; 

GETHtL3U.il>: 

GETFIELU J«l, 4 7); 

GETFItLJIJ*1» J5>; 

nt?u1Cc. POINTF R( I[ L, 0 ]) bY 4 " OflMPl 0300000" FOR NWCRDS WORDS; 

f’kKCT :* 0; 

.lEANUj :* k| L,n); 

U L EuL 2 THEN LIIlIT: = 25b LLSc. LIPiT;*1024; 

» LinIT*25fc CCikRESPONDS TO 22.5 LEG. A2IPUTH 

4 LI1If*H'24 C3RRF. SPONUS 10 22 DEG. TILT OR . 5*C ACCEL. 

J := l; 

CCPPENT ERROR IN LATA FIRST OCTAL DIGIT DIFEERENCE; 

IF ABSf RlL.J1 - K[L,J-1)> GTR LI1IT THEM 
IF Ef.RCT :* *+l GTR 4 THEN bEGIN 
*kIU(PKINTER.<“EKKCR IN READTAPE: TOO PANT BAD DATA POINTS. , 

- last one at ", 14, J); 

Go lRIOr; 

END ELSE EECIN „ , . . „ 

WIITEIPRINTER,'"ERROR IN READTAPt: 3A0 DATA POI>T AT: ROW ( , 
II, COLUrnt", 14, "). IN ARE RENCE TO ", F7.2. ", POINT ' 
F7.2, " OUT OF LIN E">, L. J. R[ L, J-1 ], R|L,J1>; 

«[ L. J 1 :» R1 L. J -11 ; 

END; 

PEAN(L 1 :* PLAN[L1 ♦ R[L.J|; 
iND JNTIL J :» * ♦ 1 GEJ NWORDS; 

XPEAN :* E E AN (1.1 / NWORDS; 

Jf : * lAl[L ]; 

f,LAN[L| :* IXNkAN - 2MP ) * CF; 

tnTE(MINlER t .< ,, PEAN VALUE OF CHANNEL ",Il." IS ".Hl.b, 

" CFl./StC.**2">,L,. , iEAN[L| ) , , tL5E 

WKITEtPRINTER.<"KEAN VALUE OF CHAINED ",I1," IS ".Fll.b, 

" RADIANS.">,L,PiAN|Ll>; 

IF l. EJL 4 THEN W RI TEI PR INTER. </>> 

D U .4 J L. J J IR| L.Jl-XPEAN) * CF UNTIu J;* * ♦ 1 GEQ NWORDS; 

END; 

S3 LEAVE; 

ERFOR: ERR :* TRUE; 

LtA /1:END RkAUTAPE; 


t nOC:Onbl:b 

OnC:On3 2: J 
b OOC:00 bb:3 

0flC:0O J b: 2 
nrcinnsi: 1 
noc; 0(117 : 2 
nnc:f0SL:3 
5 nnC:00} F; i 

one:nn A2:o 
one:OOA 2: 4 
One: OOAb ; 1 
one:00A9:n 
one:OOA 9:n 
nnc:nnA9:n 


5 

6 

6 

B 

b 


4 

REAETAPElHOC) 
3 


DOC:0149: l 
nnc-.Of A9:4 
nnc:m 19:1 
One: f'OAD: 4 
one: mpO: 3 
nnc: 0 nb 2 :2 
on: : mbE: i 
nnc:rnbF:l 
30C:03BF:4 

pne:nncl:3 
nf.Csfnci :3 
ni'C:0r.Db: 4 

o n C: 0 0 U A: l 
one:nncA:l 
nOC:OOuE:3 
nnc:rntn s 2 
n n c: o 0 ^ 2: j: 
nnc:nnE3:5 
nnc : mE6: l 
nriC:0nE7:2 
OfiC: Of'E9: 4 
nnC:03F8:l 
nnc:nOFbsn 

nncsniOAsl 
nrc: 01 11:0 
nnc:0111:4 
nnc : n 11 b:1 
nnc:nn R: 3 
nnc:011B:0 
0nc:nil8:4 

IS 012E LONG 
005:0330:4 


*< 5<<<<<<S<<< 
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I 


FhPCEOUk 1 AIJPhASE (K.I . S I. CD. HEX. CHANNELS. NY(,,SAH IN ); 

» 

F 

F^A- ARRAY R,I[1,1|.SI,CG(')].- 

il'TtCUi *1 fcA . CilANNLLi; REAL NYU; BOOLEAN SArilN; 

\ 

« 

* FCF EACH CHANNEL, J, THIS PROCEDURE CORRECTS THE TIHE-SER1ES BY 

Y 1JLTIPLYINC ITS FOURIER T.RANSFOR1 BY THE APPROPRIATE INVERSE FRE3. 

1 RESPONSE 0r ITS OUTPUT ULTtR AND THEN APPLYING THE REVERSE 

l TRANSPORT TO THE PRODUCT. Tit RtSJLTING REAL TIKE-SERIES IS THEN 
•> SOARED AND RETURNED THRCUCh ARRAY R . 

* THE CHANNEL «1 1 Irt HISTORY IS CONVERTED TO SURFACE t EVATIUN. 

» »MEN HE sAfIPLE * HOLD 1JDULES ARE iPERATING, SET SAHIN=TRJ£. 

V NOTE: PROCEDURE SINCOS HUST BE CALLED PRIOR To CAILINC ADJPHASE. 

» ADJPHASE REDDIKES FURLaAKD LEC LAPATIUN OP ARRAYS. 

Lu! " ' ’ ‘ .. . *.. * ... 

1 IT EGER HFPTS,J.K.N.NPTS; 


Kc.Al A,rt,L5, -D5,DS,P l, P'2, FI, SCAlE, SJS,SN,Sv, 2 ,TAU, TENPI, TcHPR. VAS,», 

P., W11, *12,1*0, WG2,2C, ZS; 

% 

LPTS:--( N 1 [HEX: o.-l ]; V =2»*.'1EA 

jCAuE:*2./NPrS**2; 

fi: = 3.l4l5 r 3iC5J5l>; FPPTS:*NPTS LIV 2; Sl,2 : * J * KT< 2. ); k :-l . * P 1 *N Yu/HPPTS ; 
rHJ:«T. tl; i(SECONDS) SCAN' INTERVAL BtHEEH CHANNEL SAMPLES. 

H:®!*' 1 * - A(HEhTZ) CUT-IFF f R EC 1 . ur CHANNEL fll FILTER. 

I- 2 : «= 2 f'. n ; Y i rl ERT Z ) CUT-OFF rRt*,. FOR CHANNELS 2,3,4 FILTERS. 

'* I i: * l. / (t. * PI * F1 >; WI2s»l./(2.*Pl*F2); 

l 

for J:* I STEP 1 until 4 TO IF CHANNELS. {4-J: 1 ] NEi n THFN BEGIN 
r PT(R[ J,* ] , I[ J, * |.SI,C0,1EX); El TS V2 ( ® [ J, * J, I [ J, * ', HEX > ; 

VAJ.t. :»().; 

: J TI n 0 VA.7N:»* *•< R| J,:,'| * *2 f l(J, i j**>)'SCA lE JNTi„ N:**»l GE3 HFPTS, 

IF J SOL X THEN 

Rki rs( rrintei,< #, vanian c t jp input ri.it history chan. «mi,« is «,fi2.5, 

” CF!.**2/SEC.* A 4">,J,VAkf!) ELSE 

ARlTilPRlNfil.t-VAKUNCE OF xNPJT Tilt HISTORY CHAN. 4”, 11, * IS ",E12.5, 
" RADIANS*‘2">.J,VARN); 
i 

DJ:*(J-1)*IAJ*R; SDS:=SIN(DS); CDS:=COS(DS); 

T HR j NPi'S DO BEGIN 

IF N IF.C HFpTS THEN K:«-N ELSE N:«NPTS - N; 

4 WHERE S*<I-1)*TAU*W*K; K*l,2.. 
IF N NE 5 0 AND J Ntg 1 THEN BEGIN ZC:*CS; ZS:*3N; 

If II EUL HFPTSil THEN ZS:*-SJ.'; 

C3:*ZC*ZdS*’ZS*SdS; S'): *Z3*CDS-ZC*SDS; 

ENL ELSE btGIN 


CS:*1.; SN:*f'.; END; 

SN*SIN( S)*3IM S*3S) , C5*CGS<SI * COS I S»DS). 
TE1PR: * R( J, :l] ; TE1P 1 s * 11 J, N ] ; 


O'.3:01 j 1: , 
005; or j{l: 4 
11 b: 11 j 1 ; 1 

nos:or jO: 4 
005:0 i j 0 :4 
005: 0030;., 
0P5:Of 30;4 
005 : 013 1:4 
005:0130;. 
00s.-Of 30:4 
005:01 311: 4 
005: 01 30 :1 
005:1131: 1 
005:Or jO :4 
005:01jO : 4 
01 c '30:4 

005:0130:1 
005: Of. jf : 4 
005: 11 O 3 O: , 
005:00jP;., 
005:01j1:; 
005;(r jO ;4 

3 ;H j : ni 3 n : a 

ADJPHASE IS SECKEM rOOOD 
nr>:onn : i 
noL:1000:0 

013:0,110 .1 

or o:noon: n 
10 D: HI 1 : j 
onr:TOO 3 :5 
013:0110; 1 
O n L:OOOf : j 
OOD: non : f 

000:0012:3 
0 .3:011a:, 

Of D:OPlA : 4 
010:011D:5 
Of 3: 0026:3 
113:0027:1 
100:0021 : 1 
003 : 0 I 2 r: 3 
000: 00 32:1 
010:0110:1 
000:0042:3 
003; 0 f sO: s 
i)0D:0l3l: 1 
000:0055:5 
003:1150:3 

5 000:0054:4 
OOP:005F: 1 
0iD:0i3F:1 

6 OOD:01fc2 : 5 
003:0055:3 
oot:f OCA: 1 

b OOD:0 f bA : * 


6 


1 0 D: 0 3 a D: i 
013: 005 0: . 


49 


* UNCTIONS A AND ri ARE REAL AND I MAS 1N ART PARTS OF INVERSE FILTER RESP. 

IK J El L 1 THEN BEGIN 

W02: *1 <iO:*W*K*»Il )**2 - 1.: 

IF SARIN THEN oESIN A:**0?; H:*SJ<*rfO: END 
ELSE BEGIN 

A *N32*CS - SQ2‘»0*SI; 

B:*b02*SN ♦ Sd2‘«0»CS; lND; 

IK N KQL HFI'TS THEN 

R[ 7, N): *R* TiNPH ♦ £»* T r?IP I: 

I [J,H)s»A*TEPPI * E'TEflKk; 

ENL ELSE BEGIN 

rf3:*N«K*rfI2; 

IF SAIIIN THFN b E G1N A:»l.; B:*»0; END ELSE BEGIN 

A: * C S ♦ NO*SN; B:**'C*CS - SN; END; 

IK N EJL HKPTS THtN b:*0.; 

K [ Ji N ]: = A *TEf.PR - E*TENFI ; 

I[ J.N] :*A* TE1PI ♦ d*TL1PR; 

LNL; N:* *♦I; 

end; 

• 

U.U 4 CONVERT ACCEL--RATION TRANS KO R1 TO SURKACt ELEVATION. 

IF J F|.L 1 TFiEN tELIN 

THIJ INPTS-1) DO rtEoIN 

IF N LEG HrPTS THEN 

«Os*-I./(N*v. 1**2 ELSc VO:*-I./K NPTS-M*W)**i; 
r [ l, NI: * 4 0 * A [ 1. NI; i[ 1, N): * ROM [ 1, N); 

N : * * ♦ I J 

t.nl; 

END else begin 

0IIRV2I K|J. ‘I'll 1 . * 1. flr.X ) ; FKT HI A [ J,» |, 1[ J. » |. S 1, CO. -.EX) ; 

t 

|.;*0; DO F[J,i> ) s * K [ J » N ]/NPTS UNTIL »l s * * ♦ 1 GE U NPTS; 

REPLACE POlNTcRII[J,*1) BV 0. FOR NPTS RORDS; 

INC; 

wNO; 

INI fLJPHASE; 

', » » M M U M » M M M * * . » M M M » » » M » M 


6 

7 

7 

7 


b 

7 

7 


b 

5 


5 

6 


b 

5 


5 

4 

ADJPHASc. ( 00 j ) 


000:0070:5 
OOD:f 070 : 5 
OOP:0072:0 
0 0 [): 0 7 7 o : 4 
Of D-.f r79: 4 
00o:00/A: 1 
OHD: f 07Ii: 0 
OOD:P07F:5 
00J:OOSI :5 
000:0035:3 
OOD:P0b9:i 
003:0033:4 
001:00bb:J 
OOD:OPbF:5 
OOD:0033:3 
001:0055:3 
003:0039:1 
OOD:009C:5 
OODjOOJE:1 
0(C:005F:2 
n 0 3 ; 0 0 3 K: 2 

OOD:OrA0;0 

010:0041:1 
ODD:00A5:4 
0 7 3:0 0 A o : 3 
ood:OOAE:0 
000:0733:2 
OOD:Or P4;< 
OOU:70b5:5 
ni n: robb: 2 

OOD: 00 3o:2 
1 ) 0 3: 003F : 0 

OOD:OObF:0 
i) 0 u; 01C 5:4 
OOD:00C6:1 
nnj ; mc8:i 
OOD:OOCA: 3 
IS 00D1 LOT 
005:0030:4 



50 


LAST:* 4"C5“; 

«' 13 T: * I-Ej", 

TfetbUEl* J.IoiW :* n i 1(44:1) i b t J J : 31 V 11 3 0 s 1 ]; t 7-TRX, 200 BH 
1AFEbUFj * j.IGPAaK : = n i, 1(7:1); * SUPPRESS PARITY ERROR RETRY. 

ester xnyujisi * 2./nwords; 

OPEC :* 1./I2. * FI * F STEP»; 

SISCOS! SI , CO, 1E.O; \ SET UP SIN-COS TARLE FOR FFTF 
J :« 1; 

1C PECIN 

) Rti(J) : « KSTt P * J; 3 

PEkJJ) : * 1./eK tO[J ); 

UPEo4[J) :* I orEC/J) ** 4; 

<r/0( J ):«( <3NE( J] :«.0402840996*1 FSTEP*J»»*2 >**2; 

JACOB IJ 1 :«480200./( 2.* PI»FSTEP*J)**3; 
tSO UNTIL J :* * * 1 CEO Hi 

PSI f 1 :* 0r.EC4( n ) : * ERED( 0) :« PER| 0 ) :« 0 .; 3 

XX : * (1; 

Ft ADI CARr.<2(I5,X5),E6.1,Xl,Al,X5,F6.1>,IAV.SAHIN,PMCDtV,t0RN,TLin>; 

K£AD( CASU,<b(L5.X5)>, 7A3LE1 ,Ta6lE2.TA3LE3 . PUNCHXS, PJNCH?S,PUNSHDS) ; 

It TLIP LEQ 0. THU. TLir.:*NMl/XKYtUISr; 

IE LAV HEN .RUE! P Rl N i t R, <" SPEC TR AL BANDS AVERAGED 04 LuCRUHHIC ", 

"BASIS.">) ELSE BEGIN 

Fr<EtDu)1:*2*NREC*AVEF j; 3 

RRITt(PRINTER,<12," ADJACENT FREOS. AVERAGED YIELDING ",I3, 

M DECREES OF ESlECOM It: SPECTRAL E STI PATES . "> , 

AVEFQ,FRELLOP); END; 

IF HEX LEj i AND LAV IHEN BEGIN 3 

t*RlTE(FAJNTEK,<"TIHE HISTORY TOO SHORT TO IPPLI.1ENT L0CRITHP1C" 3 
," AVERAC1N3."*); CO 10 FIMI; END; 

«RIT.I PRINTER, <"SATPLE V HOLD NODULES IN*".Lb>,SAHIN); 3 

VFITE(PRINTER,<"PACNETIC CCIKPASS IEVI AT I ON*",F6.1,XI,Ml, 

7.ACUtV,t(lhR); 

1ACDcV:* TACDEV * PI/1BD,; 

ERITH FA IATFR,<t( L5,X5)>, TABLE 1.TABLE2, TABLET', PUNCHXS.PU NCHPS.PUSCHDS); 

E 

lAPESTARlUlOYTAPL, TAPEE'IF , NNURES+l. AF ILL, FiRSTREC, EAR); 

$ 

RFC : * 1; 

I. EXT : HECIN READTAPE (MILE, REC*F IASTREC-1. CHANNELS); 

IF NOT ERR UtN ADJ PHASE I K. I, SI , CD, 1 EX, CHANNE LS, X NYO J IS T, S AHIN ) 3 

ELSE CO 70 AGAIN; 

NN:= NN ♦' 1; 

< 

K:*0; 

IF LORN tuL EAST THtN DO R(2,<]:»* ♦ HE.W2) ♦ PACDEV 

UNTIL K.:**»l CEJ N'^ORUS ELSE 
IF EOH t*,L RE ST THEN 00 R|2,X]:«* ♦ PLAN ( 2 ) - HACUEV 

UNTIL K:«"4i CEO NnORCS ELSE CO TO FXIT; 

K: *1; 

* 10NVERT ANCLE TIPE HISTORIES TO SLOPE TIPE HISTORIES CHANS. 3,4. 
irIRU NWORDS DO BEGIN 

5INAZ:«SIN(R|2,x)); COSAZ:«COf(R|2,K 1 ); 4 

TANROLL:*TAN(R (3 ,K ]); ThNHTCH:«TAN (R[4,K)); 

13SP ITCH; *C JS( R( 4, K )); 

11 3, K ): * • COSAZ*TANROLL/COSPITC(t - SIN AZ* T ANPITCH ; 

< ( 4, K ):* »SI4AZ*TANR7LL/COSPITCH - COSAZ*TANPITCH; 


nnb : m jO:» 
nP5:0fjO:4 
005:0031:J 
00 5:00 j 2:2 
005:00 35: , 
005:0037:0 
or b: 00 3A: 1 
0 0 3 :0 0 3 F: i 
003:0011:1 
005:0041:5 
035:0041:5 
003:0043:4 
005: PO<tb: 3 

005:00,3:0 
005 : 004E:0 
003:0734:3 
005:0056:i 
005:0039:5 
005:r05A: j 
015:0070:5 
005:0069: 2 
Ol5:013d:5 
005:CUBE:3 
035:0112:4 
003:01)4:1 
005:1056:3 
OP5: f 057 : Z 
005:0144:4 
003:00ftfe:3 
015 :Onn d:Z 
0O5-.O0AC: 3 
Dft5: f'f'hA: 4 
O05;00bD:2 
OOb.-OlCrt: . 
005:(OCF:, 
ori5:OU7:4 
Of5:0fE7:4 
005: 00 tb:4 
005:0 IEB: i 
005:00 tC:2 
015:00tF:2 
005:00 F4:1 
005:00)4:5 
005; OOFo : 1 
00 5:0 0 F 6:1 
(105: OOF5 :5 
005:rOFA:2 
003:00 r'J:2 
005:0101:2 
003:0114:2 
005:0105:0 
005:0115:1 
005:0105:1 
005: P10D:1 
003:0111:3 
033:0113:4 
003:0117:4 



If tORN B.JL EAST TH-N RI2.K): 1 * - f1EAN[21 - HACDE/ ELSE 

R[2,K J:«* - f!LAN{ 2 j ♦ MACDE V; 

K:“*♦ 1; 

END; 

E 4 

» 

FOR NCHANNlL:* 1 STEP 1 UNTIL 4 DO If CHANNELS.[4-NCHANNEL: 1 | NEQ 0 
ThEN DFCIN If NCHANNEL CEO 2 THEN DECIf! 

f f T ( R[NCHANNEL,*), I (NCHANNEL, * 1, SI, CO, HEX »; b 

1) ITRV2 (k| NCHANNEL, *), 11 NCHANNEL, *|. HEX); 

END; 

E CALCULATE POWER SPECTRA; CALCULATE VARIANCES. 

VP.PN: *0.; K: 1 0; 

0 J 3 r.CIN 

X :* (R|NCHANNEL, K| ** 2 ♦ I[NCHANNEL, K| ** 2>*SCALE; b 

If NCHANNEL LIE ) AND A LSS INTEGER!1./(ILIP*fSTEf)) THEN 
X: * 0.; 

A |NCHANNEL.K | : *R (NCHANNEL, K); b( RChAN'NFL, K | : *11NCHANNEL,K|; 

2 AKH :* * ♦ x; 

PCAIHCHANtEL, K ] :«= If NN EyL 1 THEN X ELSE 

psa(n:ha)nel,k| ♦ x; enj until k:=**i ceu n; 

t - 1 

;»se NlHAnmei.- 1 of atciv 

WRITE! PR INT f.R, < / < "SURfACE ELEVATION VARIANCE 1 ",Eli.b," CM * *2 "> , VA AN ); b 
WRIIEf PRINTERS "HExDINC V Ah 1MNCE 1 ",El2.b," RAD IA NS “ 2" >, V ARN ). 

Wf.ITtlRR INTER,< “N-S SLOPE V AA IANCE 1 " , E 1 2 . b ," SLUPE**2 (NO UMTS)">, 

VAIN); 

NF HE I PRINTER.* "t-W SLOPE V AH IANCE 1 ",E12.b," SLUPE**2 I NO UNITS)", 

////>,VARN ) ; 

EPD Of cases; 

4* J 

ENL; E f ND Of t FT LOOP. 

* 4 

V CALCULATE CROSS SPECTRA 

it;*1; DO CjN(K|: 1 **(A(1,N)*A|3,K) ♦ b(3,K|*B(l,K| )*SCALE 
UNI1L K: 1 *♦ 1 CEO N - 1; 

K:«l; DO ()SN[<| : »* ►( A | 3, K |*9( 1, K | - A(l, A|*S( 3,ICJ)*SCALE 
UNTIL Y : **♦1 CEO li - 1; 

K: 1 1; 00 CS21<|: 1 **(A(1»R)*A(4,K| ♦ 611,K|* d11,X |)*SCALE 

UNTIL X :**♦! CEO N - 1; 

K:»l; DO OSR|<|: l **IA|4,lC |*d(1,K| - h(l,K|*B|4,K|) r 3CALE 
UNTIL ►.:**♦! CEO N - 1; 

r.J»l; DO CNW[ N) : 1 * ♦ ( A | 3, K | * A ( 4 , K ] ♦ 8) 1. K) *B| J. K| I *SCALE 
UNTIL K.: 1 *♦ 1 CcO N - 1; 

K: -1 •’ LO y NW | K ): 1 * ♦ ( A [ 4, X J * b (3, X | - A | 3, K )■ B | •» ,K | ) * S ChLE 
UNTIL K:■ CEy N - 1; 

K: = 1; D3 PS| < | ; 1 *♦ ( A | 1, K ) * * 2 ♦ B|1,K|**2)*SCALE/0!1EC4|IC) 

UNTIL f :«**1 CEO N - 1; 
t 

ENL; 

4 3 

AUAlN:REl':***l; 

)t NEC LEJ NREC then CO TO NEAT; 

E 

t< I I: If NN ,TK 0 THEN BEGIN 

NCHANNEL: 1 !; 2s i l./NN; 3 


nnb: 011 o:3 

onb:Oil) :4 
0nb:0123:0 
005:0124:2 
00b:r1 zb:3 
105:0125:3 
005:0125:3 
00b:( 1*7:1 
O0b:('l2A:0 
00b:ri*E:2 

00b: f 1 JIM 

nnb: 01 3 1 :4 

OOd:! 1 31:4 
005: r 1 31: 4 
00b:013j:0 
003 : 01 j 3:n 
0 (■ *>: i* 1 j 7: 3 

00b:01 3d:3 
00b:013C:1 
005:0142:3 
005:0143 :i 
003:0147:3 
orb:014b:b 
005:0143:b 
00 5 :014E: 3 
0 Ob:01bD: 1 
00b:01bp.:4 
003:0160:1 
00b:017A:1 
0 0 b : 017 C: 0 
00b:01Bb:4 
Ola.-nUBsO 
00b:01bb:0 
005:0130:* 

nob;( 1 SL :2 

00b:f‘18l:2 

005:0133:2 

005:0136:5 
or.b: niSC: b 
005 : OUO : z 
00 b:01Ab: 2 
005:0143:3 
005: 01 A f: 1 
OO 5:01d J : * 
005:0113:4 
00b:01bD:1 
nr b : 01C3:3 

i)05: o 1C 7 :0 
OObsOlCC:1 
nnb:Oi:F:i 
00b: M Cf: 4 
005:012f:4 
OOb:01Cf :4 
003:0101:0 
00b:01L2:2 
nib:01D2:2 
005:0113:3 
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DJ 1EGIN 

IF 803 LEAN (CHANNELS.[4 -NCHANNEL:1|> THEN 
IT NChANNEL E(/L 1 THEN 00 BECIN 

PS[K 1 :** * Z; 

PS A [ 1, K | :* * * Z; tND 

UNTIL K:**♦! LEQ N-l ELSE 
00 PSA[NCHANNEL.K):«* * Z 
JNTIL F.: * * ♦ 1 CEO N-l ELSE 

HI EL3 T E L BESi5 N ARITE(PRINTER?<"UNSUCCESSFUl tape *EAD">); 


CO TO F INI.* 


ENU: 


CSW(K):** * Z; 
JNiW [ K j: * * * Z; 


K:*0; 

IF NN 7 T K 1 THEN UO BEGIN 

CbMtK ):=* * 7; QSN[K):*‘ * Z: 

,3a[K|:*‘ * Z; CNW[K]:*‘ l: 

EN'J UNTIL K:**♦. GEg N; 

wri rE( pKiNiEkl skip 11>; 

FOR NCHANNEL:=l.3,4 00 BE 11N 
V A r K: * 0 .; K: * 0» 

DO YARN:* * ► PSA [ NC HANSEL, < 1 UNTIL «:*»♦! CEQ N; 

IF NCHANNLL EgL 1 THEN _ „ ., „ .. 

WRITE! PNlNrEH,<«AVERAS^VARIANC t ^£j i *ij"* /> ^HREC?NCUANNEL,VARN) 

^ j t p ‘i P w l N r P R <" a VEHAk«E VARIANCE 3F ",12 ," RECOF’DS, CHAN.*", 

II,",E12.S»" SL0PE**2 (NO UNITS)",/>,HRtC,NCHANNEL,VARNI; 

ENL; 

t 

•Rir-'^PHlNrES.R//,"UNNOHHALIZEO POWER SPECTRA AVlKAGED OVER ",I2. 

" RECORDS.">,NKEC >; 

WRlTKPRiNTEk,<XSS."ACCEL.",X3."SURFACE".X4."COUP. . 

«IIE( , PKiNIEH. l "i”sUP"“x2" , -F#EQUE1C<". <2, "PERIOO". X3, "ONEOA-4" 
, Xj' "JhAN^L.S*:", 5Iin/>. F.k NCHANNEL ;* 1.1.2.3.4 DO 

IF CHANNELS. ( 4 -NCHANNEL:I 1 EgL 1 THEN NCHANNEL); 

IF CHANNELS.[4-NCHANNELjI| EQL 1 THEN PSA[NCHANNEu,J) 
jNTIL J:*‘»l CEJ N; 

ENL; 

'if CHANNELS L3S IS THEN CO TO FINI; 

ARVEFAGE TOGETHER RIUACENT FKEOLENCY BANOS. 

;:«l ; 

IF LAV THEN BECIN 

X:*1./FSTEP; Y:«-X, . . v. .on * 

00 BECIN AVtP FK | J 1 :* 1 FhEO (J J; PS ,<| | :* X. i'j]:,* * 

PSAI2,J | :** * X; PSRiJ.J):** X. PSA[4 . J J 
CSNI J ]:** * USN( J) :** * x; 

OSW[J):** * Y; CNW[J J:** * Y; 

EM UNTIL J ; * * ♦ 1 CEC LS; 


PSA[4 ,J j: * * * X; 
CSW( J1:•* * y; 

QNVij J j :** * Y; 


00b:0I07:7 
DOE : C1 i.7 :0 
nns:mo7:4 
nnb:C 1 19:2 

019:0103:0 

nnb:C1LC:b 
nnb:niLF:0 
00b;P1E1s 4 
nnb:01E4:l 
nnb:01E7:3 
nnb :01E9:3 
OATA IS OPFC LONG 
nnb:01EF :4 
(10b :nl f n : 1 
nnb:niFO:i 
(inb: 91FO : b 
nnb:niF2:0 
nnb:niF7: j 
nnb: 01 F 0 :0 
nnb: 01 FF :2 
nnb: Owl 4: > 
nnb: 02 nb:i 

0 lb:0*lC:4 

nob:Ow il:2 
0nj ; n2I2:O 
nnb:0 214: 2 
nnb:n21F :n 
nnb:C224:4 
nnb:0227:n 
nnb:O r 37:1 
nnb:nw J A:3 
nnb:nw3P:3 
mb ;n23B: 2 
FiOb: 023D: 1 
nnb:0249:1 
nnb:n24B:U 
nnb:32U:4 
nilb:n23l:J 

nib:02o2:2 

Od b: ni t»E: 1 
0'lb:''2oE:i 
onb:0230:2 
003:0230:9 
OOj 10i.3 2:4 
fiOb: 029C: 3 
0 n b : 0 2 3 C: 3 
00b: 0w9C:3 
nno:02)D:b 
or.b:029E:b 
00b:0 2)U:b 
00b:02SE:3 
003:02)F:2 
OPb:C2A3:0 
Ofib: P2AA: 0 

nnb:nibO:b 
nnb:n2Bb:2 
nnb:OwbB:b 
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VPLAVE1U:*!; tNLPT: *128; J:=b5; 

* 


4 005:02bt:l 

P05:02C0:3 


..IDi EP: VALA7EK j : «• 
1C ELGIN 

C:»lS 

lu begin 

f |*»* 
»J l J Js *• ' 


2 • 


» Mr )♦< |; 

PS| J♦ K |• 

♦ PSA j1*J*K 
J, J |* » P3A(3.J»< 

J|s»* • csnijo ); 

:•* ♦ IS.j J«K)i 
• CNN J*K i 
J|♦ JACCft(J«K|; 

<T»Q[ j Jl»* » <TNO|J»N|; 

tNi UNTIL »:>*«! CK V*l.»v«*<,; 


UNEG4| J ):•* ♦ OflEi4|J»*|; 


PiA 
CSN 
CSN 

CN. 

JACOB 


1. J 11 1 


PSA|2.J 

PLA 


gsi 

JSN 

w :.d 


4.J|:»* 
j)s»* ♦ gss 
J1 *■* ♦ os. 
J):»* • »NN 


PS A|2 <J*K |; 
PS4|4.J*<|; 
‘j*r |; 

J*K | J 
J*V ji 


KLNE|J|s ** ♦ KOBL|J*«]; 


As - 1. / ( »4L*2lFw*fsnPl; 


* i » - V i 


2l«l./VAimt Jj 

BEGIN 

OfEs l Hi "NATUHAL" NIUiiT-rtANUlD BUOY COUMUINATES AHt NORTH > REST. 
SUE LAST 4N0 JjN i4 ARE ION- CONVtNTl) N AL; SO, SINCE E. SL3Pi*-A.5 LOPE 
10P NOllULlZAIlON MULTIPLY ALL CROSS SPECTRA INVOLVING I*. SLOPE BV 
Y * -X UN EX X It ONLY AVERAGING IS TONE. 


MEuI J |: ** * 2; 

pJ4 • X; 

’SA 4,J j:* * * x; 

CSR J):** * Y; 

JACOB[JJ:** * Z; 

tND; 

pvtPEl(J):»1./tKto|J|; 

END UNTIL J:*‘ ♦ \MLAVtHj GSg ENLlPT; 


0H1G4 j J | s ** 
PSA 12,J | J «* 
C S VI i 1J • * • 
JSR[ J |:*» * 
► ONE|J ]:»* 


• i; PS|J ] <•• • xi 

• x; psai 1, J | j»• * X; 

x; oSN| J 1 1 «* » <; 

y; cn.j j ):* * • ymnv( j 1 * ; 

' Z; xT.o|J Ji■* • Z; 


ENDP T: -* * 2; 

It ENtPT Lt w hr.l THEN LO TO LCOPtK; 

». 

MU ELSE 00 BEGIN 
X :« 1; 


dj begin 

P«M)(J 1 s»• * MtaiJAK]; 

P 3[J J:* * * PS[J*K); 
PSA(l.J):** ♦ PSA(1,J + A]; 
PS A [ J , J ):>•* ♦ riBp.JAK); 
C S N [ J ): 1 * ♦ CSN | J ■» K ]; 

CSIk [ J j: s * ♦ CS. (J♦ K ] ; 

CNN[J j: * * ■» C N A | J ♦ t j ; 
JACOH(J):** ♦ JflCOP|JAt); 
< r 4 3 [ J J: = * *• <-fJO|J»KJ; 
EflU UNTIL <: = *♦! GEO AVEtQ; 


0PEG4(J]: = * ♦ 0HEG4(J♦« ); 

Pf A( 2, J j : = * ♦ P5A[ 2 ,J*K ); 
FSA(4,J]:** ♦ rSh[4 ,J*Kj; 
QS:i( J) :** ♦ ySH[ J♦ K ); 

OSW( J) :=* ♦ QS.[ J-»K |; 

ONA| J j :** ♦ y N w { JA|< j; 
KPNLj J|:* * ♦ KONE[J*K|; 


2:»1.2AVEt u ; 

BtLIN 

tPEQ|J]:-* * Z; 

P5 A [ l, J 1: * * * X; 
C S N j J |: * * * > : 
OSk[J j:* * * Y; 
JACOB|J): 1 * * Z; 

mo; 


X: >1. /(AV EtO* ESTEP); 

OfU.G4|J|:«* * Li 

; PS(J): r * * X; 

F J. ( 3, J J: ** * 

X; PSrt[4,J):** 

USNlu 1* • x; 

C5.(J):=- * Y; 

C N «j J ] : * * * Y; 

UNJ(J1:* * * Y; 

10 f 11 (J):* * * z; 

KT’niO ( J ] : “* * 


X: 


<; 


4 

5 


5 

5 


5 

4 


4 


4 


4 


003:0200; J 
005:02C2:0 
005:02:2:1 
005:02C2: 4 
005:0202: 1 
005:02:5:0 
005:t 2 CA:4 
005:0201;J 
005:02Lb:0 
005: 02 TD:2 
005:02.2:4 
005:021t:0 
n05:02Eu:2 
005:02tG:0 
005: (ut2:2 
005:02K2:2 
OO 5:02t B : j 
005;0208: i 
005: 02 r B: B 
005: 0,2tB: j 
OO 5:02tb:B 
0 0 5:0 2 1B . 3 
005:02FE:1 
005:0314:4 
005:0 31A:4 
005:0312:0 
005:0317:j 
005: 0,3,7: 3 
00i;03lB:O 
00 5: 0 j1L:3 
005:031U:3 
0 0 5: 0 311: f 
005: 0320i: t 
,105:0320:2 
005:0320:5 
005:0321:3 
005:0321:3 
005:0326:5 
005:0325:3 
005:0321:5 
005:033g:5 
005:r jJC: 1 
00 5: 0341 : j 
005: 0j 46 : 5 
005 :0 3 4 C:x 
005:0311:5 
005:0351:1 
003:0351:1 
005:0357:j 
005:0357 :j 
005:7 3 51': 0 
005: 03 b j: 5 
005:fJb9: 2 
005:036E:5 
005:03/4:2 
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AVE J Em MJ DIV AVEF'iJF ]:*l./FREy( J J; 
lM INTIE i:=*+AVEEl CtO l< - i; 

CALCULATE «AVE DIRFCTICN AND KPS BFAf HIDTil. 

IF lAV TrttS BEGIN VAuAVEFj:*l; ENJPT:*Sj; END 

tLaE CEC1N VA LA Vf F D: S AVEF y; ENORl: (NP.l; END; 


t 

N: * 1 ; 


(< 


'L LIK:DO BEGIN 

ir la; r ken j 

Clls *PSAl1.<1 

C 34 : * CNV* ( K 1 ; 

JAC OBF:(JACOB(K) 
IF tU NE'J 


ELSE J.-(l ♦ < DIV AVEF3; 

:J3:*PSA|3.<|; C4 4:*f3A(4.<] ; 

U l 3: *0SN t K 1; lU:(QS»(K]; 

C13:*CSN(K); Q34:*WNV|K) ; 

Fl1:*KONE(K ]; K 22: *KTb*0i K ) ; 

THEN Z ET A:(ARCTAN2I DU . -Q13) 

LLSE ZETA.-sSIGNI-jUUPI/*; 


ZEiA: 1 * * lbD./Fi; 

t ZEIA IS THE DIKtCTIDN OF PROPAJAHJN n'KT N-COORDINATES. 

» BETA la THE DIRECTION FROP «nICH THE HAVES ARE COPPING OEC. T. 

BETA|J|:»ltO. - ZEIA; 

K2 i * < K1 :*.l'4 r 2t)'.055G*( FS7EF*K>* * 2 >** 2 ; 
t 1 * ?I**2/JEE = .0412340 lbs. 


IF Cil NEQ 0. TFiEN 

4 k:(S yRT < ( C33 - C44l**2 ♦ 4.*CJ1**2>/<*2*CII) ELSE RR:*0.; 

r L l[J):( nK. 

C: * ( I. - RK ) / ( 1 . ♦ KR); IF C LSS C. THEN C:=C*.; 


4 

4 


t 

4 


C:*2.*AhCTANISJRTIC) ); 

CKt>r| J);*13°.*C/PI; 

OALCulATE FOURIER COEFFICIENTS OF DIRECTION FOR EA. FKtJ. 
DIREC'I I f NAL COEFFICIENTS ARE «iKT ANGLE ZETA NOI EETA. 

A 'I [ J ]:«:il/Pl J 

F11J):*11?/•F I*N1); HI|J]:»Ci4/<Pi*Kl ); 

m 2|J]:=(C33 - C44>/4 PI*K2 >; b2|J!s = 2.*C34 /(¥ I*K2 ) ; 


CALCULATE OPECIRAL POPFNTS ANO RATIOS OF POREHTS E1C 


NZc\ 1: *S 11 * J ACObK; 

MU :« w l j»JACOBK/Kll; 
•an :(*14*JACOBK/Kll; 
1.21’ :*C44*JACCbN/K^2; 

NOk :»l33*JAC0BK/K22; 
Mil :*C34* JkCC3K/K22; 


t THIS IS Non 

4 THIS IS NOI (05N). 

4 THIS IS Nin (JSE). 

4 THIS IS N2" (CC-EE) 

4 THIS IS «n: (CC-NN) 

» THIS IS Nil (CO-Nti 


UU ® Cl 4 * J A C ( bK/Kl 1; 
N2l;=ClJ*JACOBP/Kil; 
N22:*D34*JAC0BK/K22; 


4 (CC-SE) k NAB ( J]. 
t (CO-S.NI * NA7 ( U ]. 
4 t v,NE ) A NAb ( J J. 


IF NZEKO NEQ f. THEN AND 11J ] : = AFCTAN2( KOI , NlO > 
ELSE ANC.1I J J: «f). ; 

AND 11 J ] :* * * 1 b fl . / F I; 

IF NZERO ,1tO Th'FN 

ANS2I J|:*.‘j»AR , rA»2l 2.‘Ml , <N2'» - N n 2)) ELSF 
AIIF.2 [ J 

AN.2 (J ]: * * * .jn./Pl. 


JENQ11 : (Nil* * 2 ♦ NM“2; 

IF LEN0P1 NFC I'. THEN 

KNEH| J | : (< Nin*< N21 - .1021 ♦ 2. * N'll * Ml l > / DE 10M ELSE 


4 

3 

3 

3 


m5:0374:2 

nrs.-Dj Jb.-r 

005:037A:5 
005: 0 37a:5 
005:0 3 7 A:5 
nrs:r37n: i 
005: 0 37 F:4 
005 :03 7 F:4 
005:0 3bO :i 
015:0350:2 
015:0334:3 
005:038b:5 
.105:0 33C : 2 
005:0381:5 
005:0353:2 
003:0335:2 
00 5: 03 5C: 1 
005:033F’: 4 
005:P3BF:4 
Of 5:035F: 4 
0F5:03A2:1 
005:034b: 4 
005: 03AF : 4 
005:1 ja 7;2 
00 53AE: 1 
005:03AF:3 
005:031:5:0 
005:0335:2 
005: 0.3faD; 2 
005:0 jBL :2 
005:033D:2 
005:03C0;2 
005:0 3C7:5 
105:3300:j 
005:03D0:5 
005:0300:5 
005:0332:1 
005:0304:1 
005:0315:5 
005:0 3 l 7 : 4 
005:0305:3 
035:030b:2 
005:0 31b: 2 
005 ;13 DD:l 
005:0 jLF -.0 
015:0310:5 
005:0 3tO: 5 
005:0JE3:4 
005:03 Eb : 0 

n05:03EB:2 
003:03EC:0 
005:0511:5 
005:03:3:3 
005:r3 F b : 2 
005 :333 8:2 
005:03FA:2 
005:0303:0 


KNf,V| J 1 :«.).; 


MN 4 [ J 1: S I 2. » 4. * RNEW| j j/PI )/( 1 . - KNL«T[J1/Pi »; 

)E..1*12: = 3 J R r ( DE'IOMl ); 

IE IENCF2 CTR H. AND NNN [ J ] NEO n. THEN 

XAlFh[ J|:=NZERO/(2.*Pl*DKNCr2) - (NNN[J \*2. )/NNN[J) ELSE 
XAIPH|J|:-P.; 


Nni | J ); = NO 1; NA2|J|:*NlO; NA3[ J 1 : *N2*'; NA4[J):=NP2; 

4Ab( J ):*Ul; .N A b j J j :=N12; HA7[ J):«S21; NAS [ J ):*N22; 

IE 'IZEkO ell n. THEN d EG 1N 

b A l | J ): * 0.; N A 21 J ): *1 . ; SAJ[J):«1.; V.A41 J ; NAb[J]:*0.; 4 

N A b | J ) : = 1. ; Nfi 7 1 J 4A8|J|:='l.; 

ENL) ELS1 BEC1N 

NJit: T.IE W A r 10 01 THE AVERAGE HOilENTS ARE 'OMFJTED. NOT THE 4 

AVERAGE. RATIO. 

N A 1 [ J ):*’/). ZLRO; NAv[J | :**/NZENO; NA3|J): = */NZER0; 

NA4|J]:s*/NZEk0; N A b| J j : * * /NZ ERC • 

N A b [ J | :»*/!« ZERO; NA 71J ): * */4ZtR0; SAB | J ):**/NZ£RO; 

IN V l 


t..J UNTIL <:** » V A L A V l 'J CE» 14 OPT; 4 

3 

Ir (LAv AND EhDPT ELL 6b) TUN oEOlN 

K;*ENLFT; VALAVEE*:** * 2; ENLPT:*12b; END ELSE 3 

I) LAV THEN brCIN K:=ENDPT* 1; VALAVEEf:** * 2! EHIPT:* 2 ; E N D; 3 

II (LAV Ah I ENlPT LE.y NFI) THEN CO TO E.TD1R; 3 

A: * i; * E S T: * 1; 

WHILE A V E P E .* | K ] LTR .be'- TO bELIN WELT:**+f VEE * ; K:«**l; END; 

A .SiJb SEC. E L A 1C D WAVE HAS A WAVELENGTH Ly L BODY DlA.*a fT. 3 

WtST:*Nm; 4 TO PRINT LUT ALL VALLE UP U NYyUlST. 

IE' TAbLEl THEN WR 1TE ( P AI NT Er | 3K1P lj); 

IE TA1LE2 THEN BtCIN 

WRITE! FRINI Eh,<"NORMALIZEL ANL AVERAOEl CPOSS S PECTRA" , / > ); 3 

kkITE'I PRINTER,Eb); HR IT E ( PR INTFR, F in ); HR ITE ( PR INTER, H 1 ); 

WRITE!PRINTER.Pib); 

K: * 1; 

IE LAV THEN 3E3IN VALAVEE U:*1; ENLPT:*6b; 

TAE2:WR1TE( ERINTER.E'12. CO | ?.»pI *FAEy|K|,AVEPER|K |, 4 

Pj [< ]» PGA| I. < !. i'SA | 3, A |,F0.1|4,K),C3;J|k),.SN[K|,CS*(K),JSH|K), 

CN« I K 1, C'N* [ E' ). DE TA j K ), CR.eSTI K I I UNTIL K:=*+V*LAVEEy CEO 
t'JLFT); 

II LNL r T tvl L'> Tilth PtolN N:=6b; VhLA VEE y:* * * 2; tNDPT:* 12 B; 

END ELSE BEGIN K:*ENLEH-1; VALA'VEEO:** * ti ENDPT:»* * 2; END; b 

IE E.NLET LI v I..TI THEN CG TO T A h 2 ; b 

EM) lL j E 

WRITE!EK1NTER.E1/. I 0 | 2 .* PI *EPEy |Y ]. AVEFER(i ♦ IK LIV AVEEyl), 4 

ES|K|,PSF|1»K),PSA| 3.K).PSA[4.K|,CSN(K|.USN|K|,CS»|K).CSW|K]. 

CM* | F |. y N w | K |. B E T A [ 1 ♦ IK MV AVEE3)).CREST|1♦(K LIV AVEE&>|) 

UNTIL <: * * t A V E E y CE.O WiST); 

ENO; 

IE ((1ABLE2 AND TABLEJ) OR ( TABLE! AND TABLE))) THEN 3 

• Rl l'E! PR 1NTE.<| SUP ill; 


0Db:141 i :5 

nnb:0402:3 
IHa:0 4 12: J 
nn i j:040C:4 
niVa:043C:4 
nr 3: 04 PR : 1 
nnb:f 4 1 n: n 

PP b: P4 17 : 2 
n n b : n 4 1 h : o 
POD:04 lb : P 
nnb:n4lE:2 
nn b:04 t J: 4 
(103: 042 3:4 
0Pb:P4v4:b 
!IPa:n42A:4 
I'L b : P'> tE . 1 
OPb:i «cE:4 
mb:04t l .:4 
PO j : 04 2E • 4 
PO j: Pa J<<: 1 
no>:0437:b 
00 5: 04 3L : 2 
nn D :P4 dU ;i 

POb:r< )F:b 

OPb:D 4J F: b 
OPb: 04 4 1 :4 
00a: 04 4b:0 
PPL: 04 4A: 4 

POL : fl4 4C: 3 
PPj:744C:3 
nr.b: Pa sL: b 
POL:04 b4:0 
fiO b: 04 bt: n 
00b:04bb:3 
003:04 bb:f 
0Pa:14ad:P 

or b: 0.4 bb: b 
OPb: 04 bl: t. 
000:0471:3 
OP j;P4 77:2 
0Pb:047b:O 
OPb:f 47A:2 
PPb:743B:b 
00^:04 bE : l 
OPb:P 4 A A:1 
00b: 14 AF : 1 
00 3 : 04 s3: 4 
no3:04 B8 : 3 
OP 3: 04 bb : 3 
0n3:04dd;b 
nOb-04cb:3 

nn a:04 C9:4 
OP b:f 4 L E:2 
n 0 5 : 0 4 E b: 1 
00b:04 E2:1 
') n b: 0 4 r 2:1 
1)03:04:4:0 



J : = 1; <: * 1 ; 

UHILt AV EP Eh [ J ] CTR TL1M CO J := *♦ 1: 

WnILt PVEF Eh(K ] CTR . S8S LO K:***l; 

IE TAfcl.E 3 1 liEN BEGIN 

HHITEf >RINrtK, El 3, TLI'll J «R1TE( PRINTER, E 14) ; 

WHIUl PRIMER.F5); 
lr (NOT LAV) THEN 

UKITE ( ERIN1 E.R. Eb, DC ( J , AVLPER|J ) ,AO[j ) , Al[J ).A2|J ). 

o 1[ J ], 62 ( J ), BET A[J 1,CREST[J j 1 LhTIL J:=**l CEC Kl 
ELSE BEGIN 

VA LAVEE\>: * 11 ENDPT:=b5; 

UB i :*H IT E( PR INUh.Eb.DO [ J , A VEPEK) J). AO[ J ), Al (J ), A2| J ) ,B11J ).B2[J . 

BETA[J).CREST)J)] UNTIL J:=*+VALAVPEC GEu ENDPT). 
IE ENiiFT E^L 6b THEN EEC IN J:*b5; VALAVEEc:-* * 2; ENDPT:=12B; 

ENL ELSE BEGIN r , , „ . , rur .. 

J:*2NJPT»1; VALA VEr u: =* * 2; ENDPT : = * 2. END. 

IE ENCP1 LEU NP.l THEN CO TC TAB 3; END; 

LNC; 

UJR TABbE 1 -HERE IS AN UNCONDITIONAL PRINT. 

unituphinter(skip l] ); 

.KirilPuINrLR.Klb); K:=l; 

khITH ERINTEh.En.ee ( 2 . * P1 * E KEQ ( K 1 . NA 1 ( l K DIV AVEFC). 

NA2(1♦K DIV AVE)U 1» 

NPj[l»< PIV AVEE'QI.NAl I IK UIV AVEFQ ), Nm 5( 1*K D1V AVtF3), 

ANb 1 [1♦K CIV AVEEC],ANC2[1 ♦ K DIV AVEFC]. RNEN[1+K CIV AVEFQ), 
(Lu[l » t UIV AVEEC). . , 

N N N ( 1 ♦ K. UIV AVEFC) .XALPH) 1+K DIV AVEECll ENT IL N:**+nVtIC 
btt. NP.l ) ; 

i. TABIE b IS ALSO AN UNCONDITIONAL PRINT. 

*hITE(PHINlEP|SKIP 1)); . 

*«r.E< PUNIER,<"THESE MOMENT RATIOS SHOULD BE SMALL; NCHELh* . 
"NCI ♦ MO (SHOULD ECL 1.0).“,/>); 

PRINrEK.E'13); M*l; , „ . 

ARITEIEkIMEP.E lb,DO ( 2.*PI *EREQ(A ),NA6[1 +K DIV AVEFC). 

NA/ l+K UIV AVEF j j, NAB (l DIV AVEEJ), 

NA1 [ 1 ♦ K DIV AVEFC] ♦ NA2|1K DIV AVEEC)) 

DNT ID K:** + AVEEJ GEU HPl ) S 

IE PU NC H< S THtN BEGIN 
J 5 * 1 i 

a* mr bAn nts’ 

UR IT LI EUNCE-EH. < b E1« . i > , DO ( ESN [ J ), c UK ( J ). C SU [ J ), ,S» [ J ) , CN * [ J ), 
UNU( J ) ,PS[ J ), EE EQ | J j 1 UNTIL J:** + AVEE'C C£ E 
ELSE ELGIN 

VAuA/EFu:*!; t NUPT:*b5; 

FUN] :UMTl< PUNCHER, <bElDO ( CSN) J ), Q S N ( J) .CSV ( J 1.0-*» I - 1 »CNU [J ), 
j l( V ( I ], P S ( J ], E EE w I U 11 UNTIL J : - * ♦ V ALA VEE U LEU r.NDP.); 

IE ENDPT E k L Eb THEN BEGIN J:=65; VALAVEFQ:** * 2. ENDPT: = 128; 

END ELSE BELIN 

J:*ENOPr»l; VALAVEE'U:** * 2; ENDPT: 1 * * 2; END; 

. IF ENDPT LEs) N-.l THlN ,0 TO PUN1; END; 

END; iE rUNCHPS THEN BEGIN 

IF (NOT LAV) THEN BEGIN nwrl . r 

J:«l;hRlTE(PUNCHER.<bElC.2>.tO PSAjl.J] UNTIL J:**«AVEEQ oEC NM . 
J:»l;»Rirt( PUNCHER, <8EH.2>, DO PS A( 3, J) UNTIL J:***AVtFQ -EJ NMl); 
J: * 1; »'RIf F( PUNCrttH, <dEl'). 2>, DO PSA(4,U) UNTIL J: = *«-MVEFU Gr.j NM>; 

END ELSE BE)SIN FOk <:*1, J, < DO 3ECIN 

VALAVEE C:* 1» ENDPT:*65; J:*l; 

EUNiiUhllEI PUHCPIFH, *• 8H0 . * >, D6 PSA[K,J) UNT IL J :-* ♦ V ALAV E E C DEC ENDPT); 
IF ENDPT LCD 6b 1 HEN BEGIN J:*6b; V A LA V E. F J: * * * 2; lnDPT:=126; 

ENL ELSt BEGIN 

.1: -ENDPT 1; VALAVEE j:** * 2; ENDPT : * * * 2; END; 

IE r.NLE-1 LE * NP1 THEN CO 10 P U N 2 .* FND; END; 

ENL,' IE I’UNCEILS THEN BEGIN 

IE (NUT LAV) 1 HIN BEGIN 
K.-KMl UIV A V E F 0; 

J:*l ; 2RITEI PUNCHE. R,<BE10. 2>, LJ JE REjJ E»-( J-l )*RVEEJ ). SL'T<\( J ), 
LRESr|J|.A.ilJ|.Al[.J).A2(.I).Bl[Jl.b2[Jl) UNTIL J: = *M GEU M. 

END ElS e BEGIN 

VALAVCFj:•l; ENuPT: = 6 b; J: = l; 

FUI. J: U hi TL< PUNCH ER,<tEl('.2>, LO ( E RECI J 1 . b El A [ J ), 

C »I. ^ T | J )»A • (J ),Al|J].Ab|J ),Bl|J|,B2(J)l UNTIL J := * «V AL A VE FU CEC 

EriDPT); . , 

IE ENDPT EJL bb THEN BEGIN J:*bb; VALAVEP'U: = * * 2; tHJPT:=12B; 
E.»D ELSE BtGIN 

J:*LNLPT*1; VALAVE EC : s * * 2; ENDPT: = * * 2; END. 

IE ENIPT LEU KM THEN GO TO PUN 3; END; 

END; 

P IN I: ENL; < END OF PAIN bLCCM. 


ENDPT:** * 2; END; 
END; END; 


V ALAV EE'U : = * * 2; £HDPT: *1 2U : 


ENDPT:* * * 
END; 


nnb 

mia 

B neb 

6 nib 

b on 3 

4 neb 

j onb 

4 nnb 

00} 
nnb 

003 

4 nnb 

nos 

00b 

oob 

nnb 

5 onb 

b nob 

b on5 

4 nn5 

3 oob 

BLOCFinnb) IS 
STACKCODE IS SEC” 
ST ACKCODEIn19) IS 
2 on^ 

fcLOCM 002) IS 
DATA IS 


00b: 04F9: 4 c , 
005:04 EB: n 56 
005:04 EE:3 
005 : f 503:3 
nn5:0514:2 
005:0517:4 
005:051D:r 
003:05 ID: 

005:(531:5 
005:05JD:^ 
005:0543:0 
n C 5: T> S •« 4 : 3 
005:05sL:5 
005:05b9: 4 
005: 056E: l 
ni3:05jE:4 
005:0573:0 
005:(b74 : 2 
005:0574:2 
00b:0574:2 
005:0579:3 
005: 057E :4 
005:f5BE:4 
00b:05 9 4 :1 
005:(5 SC: 4 
0O5:0bA7:l 
005:05AA: 4 
n (i b: n 5 b2 :2 
005:0543:4 
005:r 5Efa: 4 
Ofo: 05 BB : 

005:(5bD:b 

00b:0bbE:t 
015:35:3:2 
n05:O5 CS:j 
91b:05D6:4 
00 5:05 IP:4 
005:0bE4:4 
0('5: 0 5 EB : 4 
0 0 5:0 5 E o : 4 
01 5:05EC:j 
M 5:15ED:1 
00b:05 EL:3 

nob: 0602 :«. 

00 5:0 oOH:J 
015:lol 1:1 
005:Ob 1i:4 
105:0027:2 
005:063b:4 
005:nbJA: 1 
005:06jA: 4 
ni5:0 o 1F:9 
005:0640:2 
0n5 : 06-.1: l 

OO5:0b4^:O 

105;06}4:1 

01j;0b3b:4 

015:0b 7 9:1 
005:OooO: / 
nOb: (‘662 : 3 
Olia: 0US4: 4 
005:0659:1 
015:0b 99:4 
O0b:ObbE:0 
005:ObA2:3 
005:0bAb:z 
OO5:06A4:1 
0 0 } : 0 o 4 5:4 
Oflb:Ob34:2 
00b:ObCA:o 
005 :0 6C A: J 
005:Ob CC:4 
005:0619:5 
005:Ob EA:3 
0O5:0bc.F:4 
0O5:0bE4:l 
O05:06P4:4 
005:('6E9:0 
OO5:06E A:2 
005: Ot P A:2 
3) IS ObEC LONG 
i SEcr.Esr com 
9) IS OObE LONG 
00^:007 7:b 
2) IS 0(94 LOliC 
fA IS 015C LONG 



57 


COMPONENTS 

The following is a list of the major components of the 
buoy, procurement source, and approximate cost in 1972. 

1. Battery: Lead-Acid 12 volts, 20 ampere hour 

capacity (6 required) 

Model GC12200 

Globe Battery Division 

P.O. Box 591 

Milwaukee, Wisconsin 53201 

2. Inverter: 12 VDC to 115 VAC 400 Hz 100 watt capacity 

Model 12400-12 (w/input-output filters) 

Nova Electric Manufacturing Co. 

263 Hillside Ave. 

Nutley, New Jersey 07100 

3. Radio Transmitter: Citizen band (Marine data channels) 

500 mw, modulated 
Model BT-109 (modified) 

Ocean Applied Research Corp. 

10475 Roselle Street 

San Diego, California 92121 

4. Visual Beacon: Flashing Xenon Lamp 

0.5 flashes/sec, self contained batteries 

Mode 1 

Ocean Applied Research Corp. 

10475 Roselle Street 


$40 ea 


$250 


$150 


$200 


San Diego, California 92121 


5. Vertical Gyro/Accelerometer Assy: Includes Honeywell $bOGO 

Model GG326C1 Quartz fiber accelerometer 
mounted on inner gimbal of Honeywell Model 
JG7044A45 vertical accelerometer 
Model GG1134AA01 

Honeywell Inc. Aerospace Division Dept 691 
2600 Ridgeway Rd. 

Minneapolis, Minnesota 55413 


6. Directional Gyroscope: North seeking 
Model DG04-0122-1 
Humphrey, Inc. 

9212 Balboa St. 

San Diego, California 

7 . Data Recorder System 

A. Data Recorder/Logic Assembly 
Model 10015 

Monitor Labs, Inc. 

10451 Roselle St. 

San Diego, California 92121 

B. Incremental Digital Tape Recorder 
Model 1387 

Precision Instrument Co. 

3 1 70 Porter Drive 
Palo Alto, California 


$700 


$4000 


$2500 


$900 


C. Low-Power Analog-to-Digital Converter 
Model ADC-12QL 
Analog Devices, Inc. 

P.O. Box 280 
Norwood, Mass. 02062 

D. Low-drift, low-power operational amplifier $50 

Model 153J 

Analog Devices, Inc. 

P.O. Box 280 
Norwood, Mass. 02062 

Buoy Hull Assy $10,000 

No Part Number 

Honeywell Marine Systems Center 
5303 Shilshole Ave., N.W. 

Seattle, Washington 98107 

Radio Transmitter Antenna: 3 needed $75 ea 

Model AW257 

Ocean Applied Research Corp. 

10475 Roselle St. 

San Diego, Ca. 92121 


